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Abstract 


In this thesis the Color Glass Condensate (CGC) framework, which describes 
quantum chromodynamics (QCD) at high energy, is applied to various scat¬ 
tering processes. Higher order corrections to the CGC evolution equations, 
known as the BK and JIMWLK equations, are also considered. 

It is shown that the leading order CGC calculations describe the experimen¬ 
tal data from electron-proton deep inelastic scattering (DIS), proton-proton 
and proton-nucleus collisions. The initial condition for the BK evolution 
equation is obtained by performing a fit to deep inelastic scattering data. 
The fit result is used as an input to calculations of single particle spectra and 
nuclear suppression in proton-proton and proton-nucleus collisions, which are 
shown to be in agreement with RHIC and LHC measurements. In particular, 
the importance of a proper description of the nuclear geometry consistently 
with the DIS data fits is emphasized, as it results in a nuclear suppression 
factor R p a which is consistent with the available experimental data. 

In addition to single particle production, the correlations between two 
hadrons at forward rapidity are computed. The RHIC measurements are 
shown to be naturally explainable in the CGC framework, and the previous 
CGC calculations are improved by including the so called inelastic and double 
parton scattering contributions. This improvement is shown to be required in 
order to get results compatible with the experimentally measured correlations. 

Exclusive vector meson production, which can be a powerful tool to study 
the gluonic structure of nuclei at small Bjorken-x, is also considered. The cross 
sections are calculated within the CGC framework in the context of a future 
electron-ion collider. In particular, the cross section for incoherent diffractive 
vector meson production is derived and a centrality estimator for this process 
is proposed. Exclusive processes are also studied in ultraperipheral heavy ion 
collisions. 
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Chapter 1 
Introduction 


Quantum Chromodynamics (QCD) is a theory describing the strong interac¬ 
tions between the constituents of hadrons, called quarks and gluons. Thanks 
to thorough experimental tests in a variety of collider experiments, QCD has 
been established as the right theory to describe the structure of matter. QCD 
has, however, turned out to be extremely difficult to solve in many situa¬ 
tions and many interesting questions related to strong interactions remain 
unanswered. 

Experimentally the QCD dynamics can be studied in many scattering 
processes. For example, the partonic structure of a proton has been studied 
precisely in electron/positron-proton collisions at the DESY-HERA accelerator 
in Germany. In current particle accelerators, such as Relativistic Heavy Ion 
Collider (RHIC) at Brookhaven, USA and the Large Hadron Collider (LHC) 
at CERN in France and Switzerland, QCD plays a dual role: it can be the 
subject of the research efforts itself or the QCD processes may be a large 
source of background reactions for other measurements (such as Higgs boson 
production) that must be subtracted. 

When colliding the heavy nuclei with each other QCD predicts [8] that a 
new state of matter called the quark-qluon plasma (QGP) is formed. The goal 
of the heavy-ion program at RHIC and at the LHC is to study the properties 
of this QCD matter in order to probe the details of QCD dynamics. These 
collisions can be described using relativistic hydrodynamics, but the initial 
condition for these simulations must be obtained using another approach. 

Atomic nuclei are more complex objects than a collection of protons and 
neutrons. For example, the partonic structure of the nucleus is not just an 
incoherent superposition of the bound nucleons. This makes the interpretation 
of the heavy ion collisions very challenging, as one has to take into account 
both initial state (the structure of the high-energy nucleus) and final state 
(the formation of the QCD matter and the interaction of particles with the 
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plasma) effects. 

In order to probe separately the initial state effects, the RHIC and the 
LHC experiments have performed proton-nucleus collisions in which one may 
not expect QGP formation. Even though these experiments were initially 
designed as reference measurements, several unexpected phenomena have been 
observed during the past couple of years. These include signs of collectivity 
for example in rapidity and azimuthal angle correlations in two-particle 
production MU- 

In this thesis an effective theory of high-energy QCD, the Color Glass 
Condensate (CGC), is discussed. In Chapter [2] we summarize how high-energy 
scattering processes, where strong color fields are relevant, can be described 
in QCD using the CGC. The framework is then applied in Chapter [3] to 
describe deep inelastic lepton-hadron scattering (DIS) and non-perturbative 
input for the phenomenological applications of the framework is obtained. 
The CGC model results for exclusive vector meson production are presented 
in Chapter |4j Finally, in Chapters [5] and [6] we present calculations for single 
and double inclusive particle production in proton-proton and proton-nucleus 
collisions and give an outlook for future directions of the research in Chapter [TJ 
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Chapter 2 

The Color Glass Condensate 


2.1 QCD at high energy 

The partonic structure of the proton can be studied, for example, in lepton- 
proton scattering experiments like deep inelastic scattering (as discussed in 
more detail in Chapter [3]). The most precise measurements of the quark 
and gluon structure of the proton come from the HERA particle accelerator, 
which collided electrons and positrons with protons in 1992-2007. 

What has been seen in these experiments is that the proton looks very 
different depending on at which scale it is measured. When the proton 
structure is probed with a photon that has a long wavelength compared with 
the proton size, a charged particle with electric charge +le is seen. The 
inner structure of the proton becomes visible when the photon wavelength is 
decreased to be of the order of the proton radius. First one observes three 
valence quarks having a fractional electric charge and carrying a fraction ~ | 
of the proton longitudinal momentum, viewed in the frame where the proton 
energy is very large. 

When the wavelength of the photon is decreased more (or the virtuality 
Q 2 is increased), a richer structure becomes visible. The photon starts to 
see a large number of sea quarks and antiquarks that carry a small fraction 
of the proton longitudinal momentum, denoted by Bjorken-rc. These quarks 
originate from gluon splittings to quark-antiquark pairs, which makes it 
possible to also determine the distribution of electrically neutral gluons that 
can not be directly probed with a photon. When the proton structure is 
measured at smaller and smaller x, more and more sea quarks (and thus, 
gluons) are seen. The extracted quark and gluon densities from the HERA 
lepton-proton data |12| are shown in Fig. 

The large gluon densities at small x are expected, as the QCD splitting 
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Figure 2.1. The parton density of the proton extracted from the HERA data. 
u v and d v are u and d quark densities, S is the sea quark and g the gluon density. 
Figure from Ref. m 


functions for the emission of a soft gluon from a quark or a gluon have a 
singularity in the limit where the gluon momentum fraction vanishes. This 
evolution, in the linear regime where gluon densities are not very large, is 
driven by the Balitsky-Fadin-Kuraev-Lipatov (BKFL) equation [13, IT]. On 
the other hand, this growth towards small momentum fractions can not 
continue indefinitely without breaking the unitarity of the theory, because the 
cross section (or interaction probability) increases when the parton density 
increases and the total interaction probability is limited by unity. This limit 
is the so called Froissart bound Hg. The growth of the gluon density can 
be limited by other partonic processes. Namely, when the gluon densities 
become sufficiently large (of the order of the inverse strong coupling constant 
l/cts), in an appropriate gauge, the gluon recombination processes gg —)• g 
become important as the probability for the recombination is proportional 
to g 2 ~ a s . This phenomenon is called the saturation of gluon distribution, 
which takes place at small x. The scale at which these non-linear effects tame 
the growth of the gluon density is referred to the saturation scale Q s (x). 

As an illustration let us consider the evolution of the gluonic structure of 
the proton as shown in Fig. 2.2 where the colorful circles represent gluons. 
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Q 2 (DGLAP) 

Figure 2.2. Evolution of the gluonic structure of the proton. 


The apparent size of the gluons is set by the scale at which they are probed, 
given by 1/Q 2 , where Q 2 is the virtuality of the probe. The change of the 
gluon density as a function of the probe virtuality can be computed from 
QCD using perturbative techniques. When Q 2 is increased, more gluons are 
seen, but as their apparent size is also smaller, the proton remains dilute. The 
evolution equations describing the Q 2 evolution are known as the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations jTBHIH] - 

On the other hand, when Q 2 is kept fixed and the longitudinal momentum 
fraction of the gluon (denoted by x ) is decreased, more and more soft gluons 
are seen as discussed earlier. As the apparent size of the gluons remains 
the same, at some point the gluons start to overlap. Eventually the gluon 
fusion gg —>■ g starts to compensate the gluon splitting g —> gg, and the 
saturation regime has been reached. This evolution in x is given by the 
Balitsky-Kovchegov equation [20l [21 j discussed in Sec. 2.4.1 As we will show 


in Chapter [3j in electron-proton scattering the Bjorken-.c probed goes like 
~ 1/s, where s is the center-of-mass energy of the process squared, and the 
saturation effects should manifest themselves in scattering experiments at 
sufficiently high energy. What is argued in this thesis is that the experimental 
data from many different scattering processes is consistent with the saturation 
picture, suggesting that the energies available at HERA, RHIC and the LHC 
are large enough for the gluon saturation effects to be visible. 

When the proton is replaced by a heavy nucleus with mass number A, 
and the nucleus is accelerated to high energy, due to the Lorentz contraction 
there will be ~ A 1 / 3 overlapping nucleons. Thus the gluon densities probed 
at same momentum fraction x are considerably different between the proton 
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and the nucleus, and the nonlinear phenomena should be visible at larger 
x. In other words, the saturation scale of the nucleus is larger than that of 
the proton at the same x. In scattering experiments with nuclei saturation 
phenomena should then be visible at lower center-of-mass energies compared to 
collisions between protons, which makes proton-nucleus and electron-nucleus 
experiments very interesting for saturation physics. 

The saturation phenomena are taken into account very naturally within 
the Color Glass Condensate framework, which is the subject of this thesis 
and presented in more detail in Sec. |2.3| It should however be noted that 


it is not the only possibility to describe saturation physics. For example, in 
Ref. |22j an initial condition for hydrodynamical description of heavy ion 
collisions is calculated from perturbative QCD using collinear factorization 
and including saturation effects. Note that the different views of saturation 
are not automatically excluding each other, as CGC is an approximation 
of QCD and other QCD based calculations can also encapsulate the same 
physics. 


2.2 QCD on the light cone 


In high-energy scattering processes the particles have velocities close to the 
speed of light and travel along the positive and negative light cone axes 
defined as 


x 


± 


1 

71 


(x°±x 3 ). 


( 2 . 1 ) 


Describing the high-energy limit of scattering processes in a quantum field 
theory becomes easier when the field theory is written using the light cone 
coordinates ( x + , x~, Xt), where the inner product is u ■ v = u + v~ + u~v + — 
Ut ■ Vt- For a detailed discussion of the quantization of field theories on the 
light cone we refer the reader to Ref. [23], and a more pedagogical description 
can be found e.g. in Ref. m- 

Let us first consider a single quark propagating along the light cone with 
momentum p, color i and spin a. This quark state can be created by operating 
on the vacuum by a creation operator b\ Q (p) as 


I in) =Nb\ ta {p) |0) =A f\p,i,a). (2.2) 

The normalization factor J\f sets the correct normalization for the single 
particle states. The fermionic operators anticommute as 

(Mp), = 5 h<M (3) (p - k), (2.3) 
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and 5( 3 \p — k) — S(p + — k + )5 > ' 2 \pT — hr)- Fourier transforming the creation 
operator to transverse coordinate space we obtain 

b la('P + i x t) = j d 2 b T e~ ipT ' XT b\ a (p + , p T ), (2.4) 

which, inverted, gives 

/ r l2 

j^ye ipT ' XT bl a {p + ,x T ). (2.5) 

Note that in the mixed transverse coordinate-longitudinal momentum space 
one obtains 


{bi, a (P + ^ x T),b^ fj (k + ,y T )} = (2'7r) 2 5 ij 5 a i 3 5(p + - k + )5 {2 \x T - yr)- ( 2 . 6 ) 


To evaluate the normalization constant J\T we require that the incoming free 
particle state (2.2) is normalized to unity. Using the anticommutator (2.3) 
we obtain 


(in|in) =M^-{k,i,a\k, i, a) = jV 2 i5 (3| (0) = V 2 J Ayje“'° = 

1 (2.7) 

where L = f dx is the size of the box in the x direction and St = J d 2 xt 
is the size of the transverse space, and we have chosen A f to be real. These, 
in principle infinite, factors are an artefact from the usage of plane waves 
instead of finite-size wave packets. As we require that the single quark states 
are normalized to unity, we get 


Af = 



( 2 . 8 ) 


The correct normalization for the single particle states is essential e.g. when 
calculating the single inclusive cross section in Sec. 


5.2.1 


An important quantity in light cone perturbation theory is the light cone 
wave function which is used to expand a state as a superposition of Fock 
states. For example, when considering the photon-hadron interaction in case 
of deep inelastic scattering (see Chapter [3]) the photon state is written as 


|7*> = l7*)o + >U’^ 5 lOT)o + ■ ■ ■ 


(2.9) 


Here the subscript 0 refers to non-interacting theory states and we neglect 
states like | qqg) which would be higher order in QCD coupling a s (or QED 
coupling CKgm) • Applying the Feynman rules of the light cone perturbation 
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theory (LCPT) from Ref. [23] one can obtain the transversally polarized 
photon wave function 


v p7*-wre = 

e o' 


e f e 


#) 


iV2—, —r^Ah(er T )(^ iT i - (1 - 2)<5 s ,±i)^ 


M 


+ m f K 0 (er T )S S}±1 S s f t±1 . (2.10) 


Similarly for the longitudinal photon the wave function is 




7 -+qq _ 


-e/e 


Qz(l - z)K„(er T )5, 


( 2 . 11 ) 


Here z is the longitudinal momentum fraction of the photon carried by the 
quark and r-r is the transverse separation of the quarks. The spins of the 
quark and antiquark are denoted by s and s', e/ is the fractional charge of 

the quark, e is the polarization vector and e = Q 2 z( 1 — z) + m 2 with rrif 
being the quark mass. This wave function is needed when the total virtual 


photon-proton cross section is calculated in Sec. 3.2 when considering deep 
inelastic scattering. Detailed derivations for these results are shown e.g. in 
Refs. [2HI 


2.3 CGC as an effective field theory 

The Color Glass Condensate (CGC) is an effective held theory that describes 
the QCD at the high-energy limit where gluon densities are assumed to be 
so large that they correspond to strong classical color fields. For a review of 
the CGC, see e.g. Refs. J2EH2B]- As the color fields are strong, of the order 
of A ~ 1/g, terms of the order gA must be resummed to all orders. The 
evolution of the color fields as a function of energy (or Bjorken-a:) is obtained 
by calculating the quantum corrections to the classical fields via non-linear 
CGC evolution equations, known as the JIMWLK and BK equations, that 
include all (gA) 71 corrections. These equations are discussed in more detail 


in Sec. 2.4. When considering scattering processes, the (gA) n contributions 
to the scattering amplitude are summed using the Wilson lines as discussed 
shortly. 

In the CGC picture the most convenient degrees of freedom are not 
quarks and gluons. To understand this, let us consider a quark moving in 
positive z direction at high energy with momentum p and scattering off a 
target consisting of a CGC (that is, off a strong color field, for example a 


heavy nucleus) as shown in Fig. 2.3. The quark scatters multiple times when 
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Figure 2.3. A quark with momentum p^ scattering multiple times off the color 
field Ai\ 


propagating through the target, but at high energy the transverse position 
of the quark can be considered to be fixed. This can be seen by noticing 
that the change of the quark transverse position during the interaction is 
A xt ~ Rkx/E , where k? is the transverse momentum obtained by the quark 
during the interactions with the target, E is the energy of the quark in the 
target rest frame and R is the longitudinal size of the target. Thus, the 
change of the quark position in the transverse plane is suppressed by the large 
energy E. This is known as an eikonal approximation. 

Now assuming that the quark momentum p has very large plus component, 
the quark-gluon vertex is proportional to gu(p) , y fJ, u(p — q) ~ gu{p)^^u{p). 
This is a Lorentz 4-vector, and as the only available vector in the problem is 
p^, we must have the quark-gluon vertex proportional to p 1 ' in this limit (the 
same result can be obtained more formally by using the Gordon identity). 
When calculating the quark scattering off the strong color field A, the vertex 
is contracted with the field A G Thus, the scattering amplitude is proportional 
to p^A^. Moreover, as has only one large component, namely p + , we get 
p^A^ = p + A~, and we expect that only the minus component of the color 
field is required to describe the scattering. From now on we work in the light 
one gauge where A + = 0. 

Let us then solve A~ from the Dirac equation 


(i$ — g4.(x))ijj(x) = 0 


( 2 . 12 ) 


using an ansatz that the solution can be written as a free particle solution 
multiplied by a factor (N c x N c matrix) V ( x ) as 


i^(x) = V(x)e w ' x u{p). 


(2.13) 


Substituting this ansatz into the equation (2.12) we obtain 


7 M [id^V(x) + p^V(x) — gA fl (x)V (#)] u(p) = 0. 


(2.14) 
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The term in the brackets must vanish with all /j. To see this, note that if 
7 ^ = 0 , then we can write 7 ^ 7 ^ = 0 , and using the anticommutator of the 
gamma matrices we get 2 a u — 7 ^ 7 ^ = 0 , where the second term vanishes 

giving a u = 0. Especially we get 

d + V{x) = -igA~(x)V{x), (2.15) 

where we used the fact that d~ = d + and that p~ ~ 0 . 

The solution to this differential equation is an exponential function. As 
V (x) is a matrix, we notice that the equation is solved with a path ordered 
expansion 


V(;x + ,x ,XT) — l—ig / d z + A ( z + ,x ,xt) 

J — OO 

( -ig ) 2 ^ 


+ 


2 ! 


dx; + / dz' + P [A ( z l+ ,x ,xt)A ( z + ,x ,xt)]+... 


(2.16) 

The path ordering P is defined such that the non-commuting fields A~ (that 
are N c x N c matrices) are arranged according to their x + component. In 


order to see more clearly that (2.16) solves Eq. (2.15) we explicitly calculate 
the second order term: 


-ig) 


2 /*aT t 


2 ! 


dz + / dz /+ P[A (z ,+ ,x ,xt)A (z + ,x ,xt)] 


-igf 


+ 


2 

{-ig) 


d^ + / d z' + A (z + ,x ,xt)A {z' + ,x ,x t ) 


J — OO 
2 px 


—OO 

+ pz f+ 

/+ 


dz' 


d z + A (z + ,x ,xt)A (z + ,x ,xt) 


' —OO 
/»*+ 


= {—ig) 2 / d z + / d z' + A {z + ,x ,x T )A {z' + ,x ,x T ), (2.17) 

J — 00 J —00 

where we have changed variables z z! in the second integral. Differentiating 


Eq. (2.17) with respect to x + gives the previous term in Eq. (2.16) multiplied 


by —igA (x). This procedure generalizes to higher order terms, and the 


solution (2.16) becomes 


V(x + ,x , xt) = P exp 


-19 


d z + A {z + ,x ,xt) 


(2.18) 
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Consider then a scattering problem where the incoming quark comes from 
negative infinity and is measured far away from the target where x + = oo. 


From Eq. (2.13) we observe that as the particle propagates through the target, 


it acquires a phase 


V (xt) = P exp 




dz + A ( z + ,x =0,x T ) 


(2.19) 


called Wilson line. Here we have set x~ = 0, as at high energy the particle 
approximatively propagates along the x + light cone, or equivalently, the wave 
function oscillates as e zp+x , and as p + is large, the wave function averages to 
zero if x~ is not approximately zero. For a similar discussion in case of QED, 
see Ref. |29| . Note that the Wilson line sums all powers of gA~ , which can 
be interpreted as allowing the quark to scatter any number of times off the 
target, and the Wilson line resums these multiple scatterings. 

Let us then consider a process where a color neutral quark-antiquark 
dipole (with quark color i and transverse coordinates Xt and y t) scatters off 
the target (strong color field). The incoming state averaged over colors is 

l in ) = (2.20) 

and when the quarks have picked up the phase factors while propagating 
through the target, the outgoing state is 

lout) = j^-VijVhlqjMqyiyT))- (2.21) 

Note that we are implicitly summing over repeated indices. Let us then 
compute the forward elastic scattering amplitude by counting the number of 
color neutral dipoles in the outgoing state: 

s = (qk(x T )qk(yT) lout) = jjrVijix^Vj'iiy^SjkSyk = -TiV(x T )V j (y T ). 

( 2 . 22 ) 

Note that, using the optical theorem, the forward elastic scattering matrix S 
can be used to calculate the total cross section. As the scattering matrix also 
contains the situation when nothing happens, it is useful to define the dipole 
amplitude 

N qq = 1 - ^ Tr V{x T )V\y T ), (2.23) 

which now includes all information about the interactions with the target. 
In practice the color field of the target off which the dipole is scattering is 
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not known, so the trace of the two Wilson lines must be averaged over the 
possible color field configurations: 


N(x T ,y T ) = 1 - -^r{TrV(x T )V*{y T )). 


(2.24) 


Note that the Wilson lines, and thus the dipole amplitude N, implicitly 
depend on the target color field and on the kinematics of the scattering 
process, especially on at which x the target is probed. As it is shown 
later in this thesis, many cross sections can be expressed in terms of the 
dipole amplitude. Thus, the relevant degrees of freedom in the Color Glass 
Condensate picture are actually the Wilson lines. In the limit of very large 
nucleus, one can derive the so called McLerran-Venugopalan model for the 
dipole amplitude [30J 



(2.25) 


where Q St0 parametrizes the characteristic transverse momentum scale of the 
gluons in the nucleus. 

To demonstrate the convenience of the Wilson lines as degrees of freedom 
let us briefly discuss the gluon distribution of a large nucleus. In collinear 
factorization it is assumed that the partons in the proton (or nucleus) carry 
zero transverse momentum. This approximation is not exactly valid at small- 
x, because the small-x partons are created in emissions of gluons with smaller 
and smaller longitudinal momentum fraction, a sequence of processes that 
can create a finite amount of transverse momentum. At small enough x the 
transverse momenta of the fluctuations can eventually be of the same order 
as the momenta of the produced hadrons in the scattering processes, see 
especially discussion in Chapter [6] where this phenomenon manifests itself in 
the production of two semihard hadrons. 

The transverse momentum dependence of the nucleus (of nucleon) gluon 
distribution can be calculated by evaluating the expectation value of the gluon 
number density operator a^(kx, k + )a(krp, k + ) in the nuclear state. This gives 
the so called Weizsacker-Williams distribution cj ) ww , which is the uninte¬ 
grated gluon distribution function of the nucleus |3T] (for a more pedagogical 
discussion, see Ref. H). It is directly related to the collinear factorization 
gluon distribution function xg via 



(2.26) 


and can be computed from the dipole amplitude N\ 



(2.27) 
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Figure 2.4. Distribution of gluons in transverse momentum space computed 
from the MV model (2.25) in arbitrary units. Most of the gluons have transverse 
momenta of the order of Q s , which is shown as a dashed line. 


Here St is the transverse area of the nucleus, and we have neglected the 
impact parameter dependence of the dipole amplitude. 

The second gluon distribution function is called the dipole gluon distribu¬ 
tion <y3 dipole , and is related to the Fourier transform of the dipole amplitude 
as 




dipole 0r, k T ) = S T 


d 2 r- 


2n 2 a s J (27t) 


T c —ik T -r T 


S(x, r T ). 


(2.28) 


The dipole gluon distribution has no number density interpretation, as it 
contains both initial and final state interactions |3l| . As we will discuss later, 
particle production cross sections in hadronic collisions and in lepton-proton 
deep inelastic scattering (DIS) are proportional to this distribution. On the 
other hand, the WW gluon distribution can be probed for example in dijet 
production that we discuss in Chapter [bj see also discussion in Ref. ITT . 

The Weizsacker-Williams unintegrated gluon distribution computed from 
the MV model and multiplied by the two-dimensional phase space factor hr, 
which gives the number of gluons with momentum hr- is shown in Fig. 2.4 


At large transverse momenta the number of gluon density drops like 1 /hr. 
At low momenta the distribution behaves like hr hi(Q s /kT), and most of the 
gluons have transverse momenta of the order of the saturation scale Q s (see 
e.g. Ref. [24]). defined as N(r^ = 2 /Q 2 ) = 1 — e~ 1//2 . If this scale is much 
larger than Aqcd, strong coupling constant is small and perturbation theory 
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Figure 2.5. Virtual photon splitting to quark-antiquark pair which subsequently 
emits a gluon. The quark and antiquark transverse coordinates are x and y, and 
z is the gluon position. 


can be used to describe the wave function of the nucleus or of the proton. 
This definition of Q 2 S is not unique, as it could as well be defined by using a 
different constant value, or e.g. from the shape of the Weizsacker-Williams 
distribution. The corresponding high-energy behaviors, for example the fact 
that the characteristic gluon transverse momentum is proportional to the 
saturation scale, are identical. 


2.4 High energy evolution equations 


2.4.1 The BK equation 

Let us find the energy, or equivalently Bjorken-:r, dependence of the dipole 
scattering amplitude. A more detailed discussion and derivation can be found 
from my MSc thesis |'2o . 

Consider a virtual photon-target scattering in the dipole frame where 
the high-energy photon fluctuates to a quark-antiquark dipole which scatters 
off the low-energy target. The forward elastic scattering amplitude for this 
process is called dipole amplitude and denoted by N. To obtain the energy 
evolution of the amplitude we consider what happens when the dipole is 
boosted to higher rapidity. The boost opens a larger phase space for the 
quark or the antiquark to emit a gluon, and after the gluon emission the full 
qqg system can interact with the target. 

The gluon emission from a quark can be computed using light cone 
perturbation theory, see discussion in Sec. 2.2 and e.g. Ref. [23] or a detailed 
calculation in my MSc thesis |25| 


The gluon can be emitted either from the 
quark or from the antiquark as shown in Fig. 2.5| Let us choose to label the 
coordinates of the quarks by Xt and y?, and the gluon transverse coordinate 
is zt- Recall that the we use the eikonal approximation and assume that the 
transverse positions are fixed during the interaction. 

The emitted gluon can be counted as being part of the dipole wave function, 
in which case we have qqg system scattering off the target. On the other 
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Figure 2.6. The emitted gluon can be seen as a part of the target wave function 
(rapidity cut at the lower dashed line) or as a part of the dipole wave function 
(rapidity cut at the upper dashedline). 


hand, the gluon can also be counted as being part of the target, in which 
case we have qq system scattering off a target with a larger gluon density, 
having a larger rapidity difference between the systems. This is illustrated in 
Fig. [2 IjJ where we also note that in addition to real correction originating 
from a gluon emission also a virtual correction coming from the wave function 
normalization requirement. The physical observables can not depend on this 
arbitrary separation of scales (whether the gluon is chosen to be a part of 
the dipole or of the target), thus equating the scattering amplitudes in both 
cases gives an evolution equation 


N q q(y + A V, r T ) = N q q(y, r T ) + -A-^Ay J d% ^ 2( ^ _ — 

x [Nqqgiy, r Ti r' T ) - N qq (y , r T )\ , (2.29) 

where N qq is the scattering amplitude for the qq dipole, and N qqg is the 
same amplitude for quark-antiquark-gluon system. The dipole amplitude is 
evaluated at rapidity y, which corresponds to Bjorken-x obtainable from the 
relation y — In 1/a;. 

Let us then work in the large-fV c limit where the emitted gluon can be seen 
as a quark-antiquark dipole. This is because the color structure of the gluon 
is approximately the same as that of the dipole, as the number of gluon color 
states is — 1 ~ iV r 2 . Thus, instead of having a qq dipole with transverse 
separation Tt — Xt — yr and a gluon at point zt, we have two gg-dipoles with 
transverse separations r' T — xt — Zt and vt — r' T = zt — Vt ■ In the mean 
field limit the expectation value for the product of two scattering matrices 
factorizes and the probability for this system not to scatter is obtained as a 
product of probabilities for the two dipoles to not to scatter: 

{S q qg(r T , r' T )) « (. S qq (r' T ))(S q q(r T - r' T )). (2.30) 


15 

















Noticing that S q q g = 1 
equation we get 


N q q g and writing the Eq. (2.29) as a differential 


d y N(r T ) = 


a s N c 

2tt 2 


d 2 r'n 


T r'^{r T — r' T Y 


x [ N (r' T ) + N ( r T - r' T ) - N{r T ) - N{r' T )N{r T - r' T )\ (2.31) 


which is the Balitsky-Kovchegov evolution equation first derived in Refs. 


12X1 • Note that we use the notation N = (N q q), see Eqs. . (2.23) and (2.24). 
The interpretation of Eq. (2.31) is clear: to get the dipole amplitude at 


higher rapidity one has to consider the scattering of two daughter dipoles 
generated by a gluon emission, and the original dipole must be removed. The 
subtraction of the non-linear term corresponds to removing double counting 
in the case where both daughter dipoles scatter off the target. 


In the small scattering amplitude limit, the non-linear term in Eq. (2.31) 


can be neglected, and one obtains the so called BFKL (Balitsky-Fadin-Kuraev- 
Lipatov) equation |13l 1T4]. whose derivation actually predates the derivation 
of the BK equation. The main feature of the BFKL evolution is that it makes 
the dipole amplitude, and the unintegrated gluon distribution function, to 
grow exponentially, which eventually violates unitarity of the theory and the 
Froissant bound. The non-linear contribution included in the BK evolution 
tames the growth at larger distances and restores unitarity. For more details, 
we refer the reader e.g. to Ref. [23J. 

The BK equation is a leading order evolution equation where the strong 
coupling constant a s is fixed. In QCD, next to leading order corrections 
are known to be significant in many processes. Thus the inclusion of NLO 
corrections to the BK equation is an important task. The full NLO BK 
equation is available (derived in Ref. |32|). and we shall discuss it more in 


Sec. 2.5 In phenomenological applications, however, the first step towards 


the full NLO evolution is to incorporate the running coupling corrections into 
the leading order BK equation and assume, that this takes into account most 
of the NLO corrections. 

The most widely used running coupling prescription is derived by Balitsky 
in Ref. |33j . In addition to that, there exist also other prescriptions to include 
the running coupling effects into the BK equation. The reason for having 
different possible running coupling corrections is that there is no unique way 
to determine which NLO terms are counted as being part of the running 
of a s , as it is a scheme-dependent choice. As an example, we mention the 
prescription derived by Kovchegov and Weigert in Ref. [34j that is shown in 
Ref. [135] to agree with the Balitsky prescription when the different scheme 
choice is taken into account. Before these perturbative calculations became 
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available, the running coupling effects were estimated by e.g. evaluating the 
strong coupling constant at the scale set by the parent dipole, see for example 
Ref. [55] . 


In phenomenological applications the Balitsky prescription is usually used, 
as it gives a slower evolution speed which is consistent with the experimental 
data 


(see discussion in Sec. 
is replaced by 


3.3). In the Balitsky prescription the BK kernel 


a s N c 


27r 2 r'^{r T — r ' T ) 2 


->■ 


a s (rl)N c 
2?r 2 


r'p {tt 

1 


rL) 2 + r! 2 


a s (r' T ) 


rpj I rp \QC S { r T Trp) 

oi s (r T - r' T ) 


(r T - r' T y 


a s (r’ T ) 


- 1 


(2.32) 


One can check that the dominant scale is the smallest of the dipole sizes \tt\, 


\r' T \ and \rx — r' T \ 


When a running coupling kernel is used one has to evaluate the strong 
coupling constant a s as a function of transverse separation. For this, we use 
an expression 

19 7T 

(2.33) 


a s (r t) 


(lliV c — 2nf) In 


4 c 2 

T IV QCD 


The uncertainty in the Fourier transform of the expression of a s from the 
momentum space to the coordinate space is parametrized by introducing 
a factor C 2 . Note that one can argue that the scale should be chosen 
as C 2 = e~ 2lE ~ 0.3152, as suggested in Ref. [51] (see also discussion in 
paper [III] and in Sec. 2.4.2). However, as we will discuss in Sec. 3.3, in order 


to get an evolution speed comparable with the experimental data C 2 must be 
adjusted properly. 

To demonstrate the effect of the BK evolution the dipole amplitude 


N(rx) is shown in Fig. 2.7 at the initial condition (MV model (2.25) with 
Q Sj o/Aqcd = 19) an( l a t higher rapidities. The Balitsky running coupling 


(2.32) is used. The general trend of the evolution can be seen from the figure: 


near the initial condition the solution approaches the asymptotic shape which 
changes very little after a few rapidity steps. Later in the evolution the 
solution propagates to smaller values of \rx\, which can be interpreted as 
having an increasing gluon density at smaller distance scales. Physically, 
this is a result of gluon splitting g —y gg. The dipole amplitude saturates 
to unity at larger dipoles, which can be seen to be a consequence of gluon 
recombination gg —>■ g processes that balances the gluon splitting. This 
interpretation in terms of gluon densities is cleaner if the BIv equation is 
written in transverse momentum space as done in Ref. (see also Ref. rat 
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Figure 2.7. Dipole amplitude from the BK equation at different rapidities using 
an MV model initial condition. 


The discussion above is valid for infinitely large targets with rotational 
invariance where the impact parameter can be neglected. In principle one 
could include the impact parameter dependence in the BK equation to describe 
e.g. the fact that if the dipole size is much larger than the impact parameter 
and the size of the target, the dipole will not scatter as the quarks are 
far away from the target. The impact parameter dependent BIv equation 
is, however, known to develop unphysical Coulomb tails that should be 
regulated by confinement scale physics (see e.g. Refs. GSHE]). As it can 
not be considered ready for phenomenological applications, only the impact 
parameter independent BK equation is used in this work. 


2.4.2 The JIMWLK equation 


The JIMWLK equation [421448] . named after Jalilian-Marian, Iancu, McLerran, 
Weigert, Leonidov and Kovner, gives the rapidity (or Bjorken-x) evolution 
of the probability distribution of the Wilson lines Wy [V ] . As the energy 
is increased, large-x color sources emit new gluons that become sources for 
further emissions and modify the probability distribution. The equation can 
be used to calculate the rapidity dependence of different correlators of Wilson 
lines. In this sense the JIMWLK equation is more general than the BK 


equation discussed in Sec. 2.4.1 which only gives the rapidity evolution of 
a dipole (correlator of two Wilson lines). Although the BK equation was 
derived before JIMWLK, it can be obtained from the JIWMLK equation in 
the mean field and large- N c limit, see for example paper mu. 
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As we will discuss in Chapters [3] and [5] for example the deep inelastic 
scattering and single inclusive particle production cross sections can be 
expressed in terms of the dipole operators only. However, when considering 
multiparticle production, one generally needs more complicated operators 
made of more than two Wilson lines. In Chapter [6]the two-particle production 
cross section is shown to require knowledge of the so called quadrupole operator 
which is a trace of four Wilson lines. The evolution for these higher-point 
functions can be obtained from the JIMWLK equation. 

The running coupling corrections to the BK equation are known to be 
very large [35] and have an important effect on phenomenological calculations. 
For the JIMWLK equation, the full NLO equation is known [49] but a 
numerical solution is still lacking. On the other hand, the leading order 
JIMWLK equation has been solved numerically (see e.g. Ref. For 

phenomenological applications it would be useful to have a way to include the 
running coupling corrections in the JIMWLK equation similarly as in the BK 
equation. What is proposed in paper mu is a running coupling prescription 
which keeps the functional form of the JIMWLK equation intact and has 
the same limiting behaviors as the BK equation with the Balitsky running 


coupling prescription shown in Eq. (2.32). 


For numerical calculations the JIMWLK equation is written as a Langevin 
equation for a single Wilson line. Following Ref. [5T] the equation can be 
written as 


dyVyM = it a Vy(x T ) 


£y(x T , Z T ) ab ' 1 £y(zT) b i + ^(xtY 


Z T 


(2.34) 


where V y is the Wilson line at rapidity y as defined in Eq. (2.19) and i — 1, 2 is 


a transverse spatial index. This corresponds to a random walk in the space of 
SU(A r c ) matrices, described by a stochastic noise term C, y and a deterministic 
part cr(xT) a ■ Averaging over the noise corresponds to averaging with the 
probability distribution Wy[K]. 

The coefficient of the stochastic term £ y is 

Eyix^Zr)^ = (^) 1/2 K(x t - z T y [1 - Ul(x T )U y (z T )] a \ (2.35) 


where U y is the Wilson line in the adjoint representation. The kernel K(xt ) 1 
is the gluon emission light cone wave function which, in the continuum limit, 
reads K(xt) 1 = Xt 1 /xt 2 ■ The noise £ is a random variable which is taken 
to be Gaussian and local in transverse coordinate and rapidity with zero 
expectation value and 

= a“W 2) (*T ~ VrWv - »')• ( 2 . 36 ) 
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The second term in the evolution equation (2.34) is a deterministic term which 
is numerically demanding to compute as it would require one to reconstruct 
the adjoint representation Wilson line. However, as discussed in paper [Illj 
and in Ref. (521. the deterministic term is needed if one wants to write the 


JIWMLK equation in the Langevin form, Eq. (2.34), as a multiplication of V 
from only one side. If the JIMLWK equation is written in a form where the 
Wilson line is multiplied from both right and left, the deterministic term is 
not needed and one obtains 


V v +dy(x T ) = exp | -i y / a ° dy j K(x t - zt) ■ (V fW ^W)} 

x Vy(x T ) exp |J K ( Xt _ Zt ) . ^( Zt )| , (2.37) 


which is equivalent to Eq. (2.34) up to order d y. In the numerical calculations 
presented in paper im equation ( |2.37[ ) (with running coupling modifications 
discussed next) is solved. 

Let us now consider the evolution of an operator consisting of Wilson 
lines. In order to calculate the operator at rapidity y + dy one has to expand 
all Wilson lines using Eq. (2.37) up to order dy, or equivalently £ 2 , and take 
the expectation value over the noise term £. Only the contribution ~ a s 
remains in the evolution equation, which physically corresponds to the fact 
that in the leading order JIMWLK equation the evolution is obtained by 
calculating an emission of a single gluon. Consider now an evolution step 
where two noise terms C,(ux) and C,(vx ) are contracted. The delta function 
in the correlator (2.36) can be now interpreted as follows. First, we note 
that if C,(ut) corresponds to the gluon emission at coordinate Ux in the 
amplitude, then £,(i’x) is the absorption of a gluon at vx in the complex 
conjugate amplitude. Let us denote the momentum of the emitted gluon by 
kx- Now, the amplitude is proportional to e lkT ' UT , and the complex conjugate 
amplitude to e~ lkT ' VT . When all possible gluon momenta are integrated over, 
the delta function 8^ 2 \ux — Vt) is obtained, which is part of the correlator 


( [2736] ). 

To include the running coupling we refer to a general result in gauge 
theories that the beta function can be computed by considering higher order 
corrections to the gluon propagator. Note that this is exactly what is done 
e.g. in Ref. [3_3j when the Balitsky running coupling prescription for the BK 
equation is derived by taking into account the quark loop corrections to the 
gluon propagator. Now the only scale available for the running a s is the 
transverse momentum of the emitted gluon kx, so we propose to replace the 
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fixed coupling correlator 


(Z(xT) a ’ i Z(yr) b ’ i ) = a s 6 ab 6^ 


d 2 h 


(2vr 


T gikT-ixT—yr) 


by 


<> n(x T r‘ri(VT ) bJ > = i“W J 


d 2 k 


(2>r) 


]Le ik T-{x T -yT) 


a s {kj 


(2.38) 


(2.39) 


It is shown in paper [ITI] that if the BK equation is derived from the JIMWLK 


equation with the correlator (2.39), one obtains an equation that has the 


same kernel as the BK equation with the Balitsky running coupling in the 
limit where the parent dipole or one of the daughter dipoles is very small. 

Let us then present our numerical results for the JIMWLK equation with 
running coupling. The equation is solved using the algorithm presented 
in Ref. |5Uj. We compare our running coupling results with the JIMWLK 
equation solved using the “square root” running coupling. The square root 
coupling is the simplest modification to the JIMWLK equation which includes 
the running of a s . In this prescription the coupling yTR is evaluated at the 
scale which is the argument of the kernel K(xr)- For comparison the BK 
equation is also solved using the Balitsky running coupling. See paper [ ITT] 
for details. 

The strong coupling constant in the transverse coordinate space is evalu¬ 
ated as 

47T 

Oi s {r T ) = -yy-i-( 2 - 40 ) 


(3 In 


a qcd 


+ 


in c 

4 e - 2 TE ^ c 


r 2 A 2 
r T ii QCD 


with c = 0.2, fi o/Aq CD = 2.5 and /3 = ylV c — |nf. The difference to the 


expression (2.33) used in the previous section is that Eq. (2.40) has a smoother 


infrared cutoff parametrized by the constant c. Also the scale at which the 
coupling is evaluated is not a fit parameter. Instead, it is given by 4e _27e /r^ 
which is taken from the explicit Fourier transform of the kernel calculated e.g. 
in Ref. [33]. As we will discuss in Chapter [3j the BK fits to experimental data 
suggest the scale to be 20/rf. for A QCD rs_/ 200 MeV, unless the evolution is 
slowed down by other effects (see discussion of the next to leading order BK 


equation in Sec. 2.5) 


Figure 2.8 shows dipole amplitude N (which is a correlator of two Wilson 


lines) obtained by solving the JIMWLK equation with the proposed noise 
running coupling, Eq. (2.39). The result is compared with the solution to 


the BK equation with the Balitsky running coupling. We find that the BK 
equation leads to a slower evolution speed and that the shapes of the solutions 
are slightly different. To better see the effect of the noise running coupling on 
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Figure 2.8. Evolution of the dipole am¬ 
plitude N from the JIMWLK equation 
with running coupling compared to the 
BK equation. The lines show the ampli¬ 
tude at intervals of 2 units in rapidity. 
Figure from paper mu. 


Figure 2.9. Evolution of the dipole 
amplitude using the proposed “noise” 
running coupling compared with the re¬ 
sults using the square root prescription. 
The lines show the amplitude at inter¬ 
vals of 4 units in rapidity. Figure from 
paper m- 


the JIMWLK equation we show in Fig. 2.9 the JIMWLK solutions obtained 
using our noise running coupling and the square root coupling prescriptions. 
Changing the running coupling prescription to the noise coupling reduces the 
evolution speed but leaves the shape of the solution roughly the same. 

To characterize the evolution speed we show in Fig. |2.10 the evolution of 
the saturation scale Q s defined as 


A = 


dlnQl 
d y 


(2.41) 


where N(r^ = 2/Q 2 ) = 1 — e 1//2 . The calculation is done with two different 


values for the infrared freezing parameter c in Eq. (2.40). The dependence on 


this constant is quite different between the BK and JIMWLK calculations at 
small Q s . The BK equation is shown to give the smallest evolution speed, and 
the Noise running coupling prescription slows down the JIMWLK evolution 
at all saturation scales. At larger Q s the evolution speed from the JIMWLK 
equation drops below that of the BK, which is an artifact from Q s becoming 
of the same order as the lattice ultraviolet cutoff. 

Finally, in order to demonstrate why the strong coupling is evaluated at 
scale set by 4 e~ 2lE jr\ we solve the JIMWLK equation using two different 
running coupling prescriptions. First is the same square root prescription 
used before, where the evolution kernels are in the form y/ajrr)K[ tt) ■ This 
is denoted by r-JoT s . We compare this with the momentum space square 
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Figure 2.10. Evolution speed of the 
saturation scale as a function of Q s cal¬ 
culated using the noise and square root 
prescriptions and compared to the BK 
results. The thin lines are obtained by 
using a smoother infrared freezeing with 
c = 1.5 instead of c = 0.2. Figure from 
paper HI. 


Figure 2.11. Evolution of the scat¬ 
tering amplitude with the square root 
running coupling evaluated in coordi¬ 
nate space at scale 4 e~ 2lE /r^ and in 
momentum space at scale k'^. Figure 
from paper III'. 


root prescription, denoted by ky/oT s . This is obtained by noticing that in 
the JIMWLK equation of a dipole one only has dot products of kernels K. 
These dot products can be Fourier transformed into the momentum space, 
and the coupling constant can then be evaluated as a function of transverse 
momentum. We write the dot product as 

V^sK(xt) -v = -ij (2.42) 


where Vr in practice is the kernel K with the same or different argument. 
When evaluating the strong coupling in momentum space 4e _2TB /(r|A.Q CD ) 
is replaced by k^/ Aq CD in Eq. (2.40). The obtained dipole amplitudes are 
shown in Fig. 2.11 and the agreement between the results obtained with 


different running couplings are remarkably similar. This justifies the usage of 
the constant 4e^ 27s when evaluating the transverse coordinate space strong 
coupling constant. 


2.5 The BK equation at next to leading order 

The BK equation describes the energy evolution of the dipole amplitude at 
leading order accuracy. As the dipole amplitude encodes the relevant degrees 
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of freedom to describe a scattering process, the BK evolution also describes 
the energy dependence of many observables such as single inclusive particle 
production and deep inelastic scattering. In paper Em we studied the BK 
equation at next to leading order (NLO) accuracy. 

The leading order calculations can give a good physical description of a 
scattering process, but the next to leading order corrections to the perturbative 
calculations of cross sections can be numerically large. Thus it is important 
to perform CGC calculations at next to leading order accuracy in a s and 
compare the results quantitatively to experimental data. For example, the 
single inclusive cross section is already known at NLO accuracy |531 - t3B] . 
and similar calculations for the deep inelastic scattering cross section also 
exist [571 158]. A crucial ingredient in phenomenological NLO calculations is 
an NLO evolved dipole amplitude. 

The next to leading order BK equation has been derived in [32|, but a 
numerical solution for the equation has not been available previously (for the 
NLO BFKL equation [o9l l6l] a solution exists [62]). The equation can be 
written for the scattering matrix S = 1 — N as 


a „S(r) = ^Ki ® [SpOStr) - S(r)] 

+ ® [S(A')S(JT - zt ') s ( y ') - s(x)S(r)] 

+ ® S(V)[S(.V) - S(A)]. (2.43) 

Here convolutions <S> are calculated by integrating over the transverse coordi¬ 
nates of the emitted gluons Zt and zt', and we use a notation r = |xt — Ut |, 
X = | xt — zt\, y — | Ut — zr\, X' = \xt — zt'\ and Y' = \yr — zt \ ■ The 
kernels are 
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(2.45) 

(2.46) 
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The kernel A\ consists of a LO BK kernel r 2 /(X 2 Y 2 ) and an NLO correction. 
Part of the NLO corrections, especially the term involving the renormalization 
scale /j 2 , should be absorbed into the running of a s . What other terms are 
absorbed into the coupling is a scheme choice, and here we adapt the choice 
derived in |33] and replace the terms including the beta function coefficient 
(3 = ylV c — I'Uf by the Balitsky running coupling given in Eq. (2.32). The 
other a 2 terms we choose to evaluate at the scale set by the parent dipole r 2 , 
as it is the only available external scale. Thus, the kernel K\ is written as 
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(2.47) 


and the coupling constants multiplying kernels K 2 and Kf are replaced by 
a s (r). For the running coupling we use the same expression as in the case of 


JIMWLK analysis in Sec. 2.4.2, Eq. (2.40). 


The equation (2.43) is written in the large- N c limit where the evolution 


equation can be written in terms of dipole operators only. If finite-W correc¬ 
tions were included, the evaluation of the derivative d y S(r ) would require one 
to evaluate correlators of up to six Wilson lines. Without a numerical solution 
to the NLO JIMWLK equation [49,, [63], which is currently not available, a 
possible way to evaluate the higher point functions would be to use the so 
called Gaussian approximation which allows one to write any higher point 
function in terms of the dipole as described in Ref. [64]. As the finite- N c 
corrections to the leading order BK equation are known to be very small [65], 
the inclusion of the N c suppressed terms to the NLO BK equation was not 
attempted in paper m- 


The Wilson lines are by definition (see Eq. (2.19)) conformally invariant, 
so the evolution equation should not break this invariance in a conformal 
field theory. In QCD the conformal invariance is broken by the running 
coupling effects, but if the running of a s is not taken into account, the NLO 
BK equation should be conformally invariant. Let us explicitly check the 
invariance. First, the Wilson lines are clearly invariant under translations 
and rotations. The only non-trivial symmetry that has to be checked is the 
invariance under the spatial inversion. Consider a transformation x M —> x^/x 2 
with respect to a point with x~ = 0. Now ( x + ,Xt ) 2 = -Xy 2 . The Wilson 
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line now transforms as 


V(xt) = P exp 
—» P exp 
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dx + / x + x T 
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x T 
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(2.48) 


which shows the invariance. 


However, the NLO BK equation (2.43) is not conformally invariant. The 


running coupling part (proportional to /3) is expected to break the invariance, 
but in addition to that the double logarithmic term In X 2 /r 2 \nY 2 /r 2 


m 


Eq. (2.43) is not invariant in conformal transformations. Before proving 


that, let us first show that the leading order BK equation is indeed invariant. 
Again, the only noil-trivial transformation is inversion. To show this, it is 
easiest to write the two-dimensional vectors as complex numbers x — X\ + ix 2 , 
y = yi + iy 2 and z = Z\ + iZ 2 ■ Now, under an inversion all vectors x,y,z 
transform as 

zi + iz 2 z 1 


z = zi + iz 2 -> 


zz 


(2.49) 
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where z = Z\ — iz 2 is the complex conjugate of 2 . Now the squared distances 
transform as 
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Using this result one can directly transform the kernel: 
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As the integration measure transforms as d 2 z d 2 z/|^| 4 , the leading order 
BK equation is found to be invariant under the inversion. 

On the other hand, the double logarithmic term transforms as 
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and is not invariant under the inversion. Similarly one can check that the 
other combinations of transverse coordinates in the kernels K 2 and I\ j are 
invariant under inversions. 

The conformal symmetry breaking is a consequence of having a cutoff in 
the longitudinal direction when deriving the NLO BK equation. This cutoff 
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violates the conformal symmetry. If the NLO BK equation is derived in the 
fully conformal N = 4 supersymmetric Yang-Mills (SYM) theory (which is 
similar to QCD but has no running coupling), the resulting equation still 
has the double logarithmic term and is not conformally invariant | 66 |. This 
confirms that the conformal breaking is an artifact of the use of a cutoff which 
does not respect the required symmetries. 

To solve this problem it was proposed in Ref. [SB] to write the equation 
in terms of the conformal dipole S(r) coni , defined as 

S(r) coni = S(r) - J dh— In — [S(X)S(Y) - S(r )]. (2.53) 


Here a is an arbitrary dimensional constant which will cancel from the 
evolution equation. When the NLO BK equation is derived for the conformal 
dipole, the double logarithmic term disappears from the equation, and an 
additional term 


2 r 2 


X 2 Y' 2 (z — zXf- ln X ,2 Y 2 ^ 2 ' 54 ' ) 

appears in the kernel K 2 - The conformal invariance is then restored if running 
of a s is not taken into account. In N = 4 SYM the corresponding equation is 
fully conformal. 

The first numerical solution to the NLO BK equation is presented in 
paper m- As an initial condition in this analysis we used a modified 


_/\2 


McLerran-Venugopalan model (2.25) with an anomalous dimension 7 


N(r, y = 0) = 1 — exp 


(r 2 T Ql 0 ) 


ln 


r A 


+ e 


QCD 


(2.55) 


where the anomalous dimension controls the power-like tail of the dipole 
amplitude for small dipoles. This parametrization is used in leading order 


(with Balitsky running coupling corrections from Eq. (2.32)) fits to DIS data 
e.g. in Ref. 153 and in paper jV]. There is, however, no particular reason 
why a full NLO fit would prefer the same values for the fit parameters. Thus 
we only explore the general properties of the evolution equation with different 
parametrizations for the initial condition. 

As we will demonstrate later, the dipole amplitude may turn negative 
during the evolution and may not satisfy the requirement N (r) —y 0 when r —y 
0, which follows from the definition (2.24). In this case also the convolution 


with Ki would not coverge. To avoid this problem we freeze N(r) — 0 in the 
region where it would become negative. 

The dipole amplitudes obtained by solving the NLO BK equation for 
non-conformal and conformal dipoles are shown in Fig. 2.12. The results are 
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Qs , o/Aqcd = 19,7 = 0-6 


Qs , o/Aqcd = 19,7 = 




Figure 2.12. Dipole amplitude and 
conformal dipole amplitude at initial 
condition and at rapidities y = 5 and 
y = 10 (from riht to left) compared to 
leading order BK equation. Figure from 
paper m 


Figure 2.13. Evolution speed of the 
saturation scale for conformal and non- 
conformal dipoles as a function of the 
saturation scale compared to the LO BK 
equation. Figure from paper m- 


compared with the leading order solution, and as an initial condition we use 
parameters Q s , 0 /Aqcd 19,7 = 0-6 chosen such that the dipole amplitude 
increases at all dipole sizes throughout the evolution considered here. We 
observe that the NLO corrections slow down the evolution speed of both 
NLO dipoles equally compared to the leading order evolution. The shape of 
the dipole amplitude, on the other hand, remains roughly unchanged. The 
leading order BK equation is solved using the Balitsky running coupling 
prescription [33] in order to get comparable results. 

To quantify more precisely the effect of the NLO corrections we study the 
evolution speed of the saturation scale A, defined as 


d In Q 2 S 
d y 


(2.56) 


where the saturation scale is defined as 


N(r 2 = 2/Ql) = 1 - e 


- 1/2 


(2.57) 


Note that Q s is different from Q S) o at the initial condition: with Q s ,o/Aqcd ~ 
19 we get Qs/Aqcd ~ 40. The results are shown in Fig. 2.13 The NLO 


evolution for conformal and non-conformal dipoles is roughly equally fast, 
and significantly slower than the evolution speed at leading order. The 
slower evolution speed at NLO can somewhat be expected from the leading 
order fits where the scale at which the transverse coordinate space strong 
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- Qs, o/Aqcd = 4 

— Qs, o/Aqcd = 19 

- Qa, o/Aqcd = 13 

. Qa, o/Aqcd = 26 

77 Jvr 


Figure 2.14. Logarithmic derivative of the dipole amplitude (evolution speed) 
at initial condition with different values for the saturation scale and the anoma¬ 
lous dimension. Figure from paper m- 


coupling constant is evaluated is fitted to the data. In order to get a slow 
enough evolution speed comparable with the experimental data this scale, 
parametrized by C 2 in Eq. (2.33), must be taken to be relatively large. See 
discussion in papers mm 

To study the evolution of the dipole amplitude as a function of the dipole 
size we show in Fig. 2.14 the logarithmic derivative of the dipole amplitude 
d y N(r)/N(r ) at the initial condition with different values for the anomalous 
dimension 7 . Generally one would expect the dipole amplitude to increase 
with rapidity, corresponding to the physical picture of having more gluons 
inside the hadron at small x ~ e~ y . What we find is that for small initial 
saturation scale Q s $ the NLO corrections are so large around r ~ 1/Q S (which 
is the dominant scale) that the solution propagates in the “wrong” direction. 

In order to interpret the logarithmic derivative shown in Fig. |2.14 note that 
if dyN/N has a constant positive value, the amplitude grows exponentially 
in rapidity and propagates to smaller dipole sizes. This is what happens 
with initial anomalous dimension 7 = 0 . 6 , and marginally also with 7 = 0 . 8 . 
With 7 = 1.0, on the other hand, we have d y N/N ~ lnr, which drives the 
amplitude towards a steeper shape. Eventually this leads to a singularity 
in the logarithmic derivative and N(r) changing its sign at finite r. Similar 
results were also found in Ref. [58j . As shown in paper m, the shape of 
dyN/N does not change significantly during the evolution. In other words, 
the evolution is sensitive to the initial condition still after a few rapidity steps. 
It is also shown that the conformal dipole has a similar evolution, but the 
amplitude crosses zero at smaller r, meaning that with the conformal dipole 
the “negativity problem” is less severe, but it still exists. 

In order to study the evolution of the dipole amplitude shape we calculate 
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Qs , o/Aqcd = 19 


Qs , o/Aqcd = 19 


Qs , o/Aqcd = 19 



Figure 2.15. Anomalous dimension 7 ( 7 -) as a function of dipole size at rapidi¬ 
ties y = 1,5 and 30 with different initial anomalous dimensions. The solid lines 
are the solutions of the LO equation and the different dashed lines correspond to 
NLO evolution of non-conformal dipole with different initial anomalous dimen¬ 
sions. Figure from paper Ref. ISE- 


the anomalous dimension as a function dipole size, defined as 


7 ( r ) = 


d In N (r) 
d In r 2 


(2.58) 


The results are shown in Fig. 2.15 where the anomalous dimension 7 (r) is 
evaluated at rapidities y — 1, y — 5 and at y = 30 with different initial 
anomalous dimensions. The fact that the solution becomes unstable is clearly 
visible with 7 = 1 where the anomalous dimension grows very rapidly already 
at y = 5. With 7 = 0.8 much longer evolution in rapidity is needed before 
the unstable region is reached. With small anomalous dimension in the 
initial condition (here 7 = 0 .6), the unstable behavior is not observed within 
the studied evolution range. Note that for 7 = 1.0 the solution does not 
evolve significantly from y = 5 to y = 30 as the evolution is dominated by 
region where the amplitude would be negative we have frozen N(r) = 0. The 
conclusions for the conformal dipole are similar. 

To study where the unstable behavior originates from we calculate con¬ 
tributions from different terms in the NLO evolution equation to the log¬ 


arithmic derivative of the dipole amplitude, d y N/N. In Fig. 2.16 we show 
separately the leading order contribution, the NLO contribution originating 
from the double logarithmic term (which breaks the conformal invariance) 
~ a 2 \nX 2 /r 2 \nY 2 /r 2 and the other NLO contributions. With large anoma¬ 
lous dimension in the initial condition (especially with 7 = 1 . 0 ) the double 
logarithmic term drives the evolution speed and is responsible for turning the 
dipole amplitude negative. With smaller initial anomalous dimension 7 = 0.6 
both NLO contributions approach zero in the small dipole limit, and 7 = 0.8 
is again a marginal case. 
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Figure 2.16. Evolution speed of the dipole amplitude at the initial condition. 
Shown are separately the full NLO and LO results and the NLO contributions 
from the conformal (no double logarithmic) and non-conformal (only double 
logarithmic) terms. Figure from paper m ■ 
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Figure 2.17. Evolution speed of the conformal dipole amplitude at the initial 
condition. Shown are separately the full NLO and LO results and the contri¬ 
butions from the In r 2 term and from the other NLO terms. Figure from pa¬ 
per |VII| . 


I 11 the evolution equation of the conformal dipole the double logarithmic 
term is absent, but a new logarithm ~ In r 2 is present. Note that the only 
other logarithm is ~ In X 2 Y' 2 /(X' 2 Y 2 ) which explicitly vanishes in the limit 
of zero parent dipole size. The NLO contributions coming from the a 2 In r 2 
term and the other ~ a 2 contributions are shown separately in Fig. 2.17 


In this case the aj In r 2 is responsible for driving the amplitude negative at 
small dipoles if the initial anomalous dimension is large. The definition of 
the conformal dipole (2.53) roughly speaking deletes the problematic double 
logarithmic term from the evolution equation and introduces a new logarithm 
that makes the evolution unstable in a same way, except that the unstable 
region is approached only at smaller dipoles. 

Based on these results and paper m we can conclude that a better 
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theoretical understanding of the NLO BK equation is needed berore it can 
e.g. be convoluted with the NLO impact factors E3EE] and applied to DIS 
processes. Possible solutions could be to include something like a kinematical 
constraint that is suggested for the leading order BK equation in e.g. Ref. |69|. 
A possible resummation procedure which takes properly into account the time 
ordering of the subsequent gluon emissions and sums the double logarithmic 
contributions to all orders has recently been proposed in Ref. [SSj. 

The strong dependence on the anomalous dimension at the initial condition 
is perhaps not surprising. As discussed e.g. in Refs. HQH72], the leading 
order BK equation (without running coupling) has a very different behavior 
depending on the steepness of the initial condition. In particular, as shown 
in Ref. |7U] . if the anomalous dimension of the initial condition in the Mellin 
space is small, then this anomalous dimension sets the evolution speed of the 
solution. This is known as the pushed front solution in the literature. On the 
other hand, if the anomalous dimension is larger than the critical value, then 
the critical anomalous dimension sets the evolution speed to the universal 
value in the so called pulled front region. 
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Chapter 3 

Deep inelastic scattering 

3.1 Probing the proton structure with leptons 


The proton structure can be studied by scattering (virtual) photons off it. As 
a photon source, a lepton (electron or positron) beam is usually used. 


Let us consider a process shown in Fig. 3.1 [ where an incoming lepton 


(momentum t) interacts with the proton by exchanging a photon with mo¬ 
mentum q. In order to see the inner structure of the proton, and not just 
one electric charge with radius ~ 1 fm, the photon virtuality Q 2 = —q 2 must 
be increased. In that case, the momentum transfer becomes very large, and 
the proton breaks up into fragments denoted by X with an invariant mass 
W 2 = (P + q) 2 . This process is called deep inelastic scattering (DIS). 


The DIS kinematics, when the target can either be a proton (A = 1) or a 
nucleus with mass number A, can be expressed using the Lorentz invariant 



Figure 3.1. Deep inelastic scattering kinematics 
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variables 


(3.1) 

(3.2) 


q 2 = ~ q 2 = 

AQ 2 
X 2 P-q 
P ■ q 


In the infinite momentum frame, where the proton momentum is very large, 
x can be interpreted as the fraction of the proton longitudinal momentum 
carried by the quark taking part in the scattering process. The virtuality 
of the photon Q 2 sets the distance scale probed in the collision, and the 
inelasticity y gives the fraction of the electron energy carried by the photon 
in the target rest frame. 

Photon emission from a lepton (which is a point-like particle) is well known 
and can be computed accurately from Quantum Electrodynamics (QED), but 
as the proton structure is a priori unknown, the photon-proton interaction 
can not be computed from first principles using perturbative QCD. However, 
one can parametrize the photon-hadron coupling by introducing the most 
general form for the hadronic tensor. Taking into account the symmetry 
requirements, one can eventually write the the invariant cross section for 
the inelastic electron-proton scattering in terms of two dimensionless proton 
structure functions F\ and F 2 , see e.g. Ref. [73]. The total cross section then 
becomes 


(3.3) 


dcr 

dxdQ 2 


4vra em 1 
Q 4 x 


xy 2 F 1 (x,Q 2 ) + 



F 2 (x,Q 2 ) , 


(3.4) 


where M is the proton mass and s = (P + £) 2 is the center-of-mass energy of 
the proton-lepton system. 

The QCD part of the scattering is now included in the structure functions 
F\ and F 2 . To leave out the well-known QED part, the experimental mea¬ 
surements are often reported for the structure functions or for the reduced 
cross section defined as 

a r (y, x,Q 2 ) = F 2 (x,Q 2 ) - ~ y)2 F L (x, Q 2 ). (3.5) 

Here the longitudinal structure function is F Jj — F 2 — 2 xF\. It measures the 
violation of the Callan-Gross relation which states that Fi = 0 if the proton 
consists of only spin-| fermions (quarks). Thus, Fl gives in principle a more 
direct measure of the gluonic structure of the proton. Unfortunately, the 
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precision of the current experimental Fl data from HERA is very limited IB 
IB compared to measurements of F 2 alone [12] . 

The structure function F 2 can be seen as a measure of parton density, as 
at leading order one can show that m 

F 2 (x,Q 2 ) = y(3.6) 

i 

where e* is the charge of quark type i and is the quark distribution function. 
Note that at leading order in perturbative QCD the structure function is 
independent of Q 2 . When higher order QCD corrections are included, the Q 2 
dependence of F 2 can be obtained by deriving the so called DGLAP evolution 
equation [I6l - fl9] . 

The journey to measure the internal structure of the proton via DIS began 
at SLAC experiment in 1960s where the partonic structure of the proton was 
seen. The theoretical foundations of the parton model were derived around 
the same time in Refs. USED- Currently, the most precise information on 
the proton structure is obtained from measurements done at DESY-HERA 
electron/positron-proton collider from 1992 to 2007. The HI and ZEUS 
experiments from HERA have published combined results for the proton 
structure functions and the for the reduced cross section [ 12 ] , 


3.2 DIS in dipole picture 

Deep Inelastic Scattering at small x can be described within the Color Glass 
Condensate framework by looking at the process in the dipole picture. There, 
in the frame where the proton is at rest the lifetime of the 7 * —y qq quantum 
fluctuations is much larger than the typical timescale of the interaction given 
that x <C 1/(■ mR ) where m and R are the target mass and radius (in the target 
rest frame), respectively. The process then looks such that first the incoming 
virtual photon splits to quark-antiquark color dipole which subsequently 
scatters off the target proton [23] . 

The total virtual photon-proton cross section is (see e.g. ra , or more 
pedagogical discussion in Ref. [23]) 

< r T,L(rT,x) = '22 f d? r T f ^[^*^} f TL (r T ,z)a^ t (r T ,x), (3.7) 

j J J 0 

where a^ t is the total dipole-proton cross section. It can be obtained from the 
dipole amplitude N which is the imaginary part of the dipole-target scattering 
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amplitude by using the optical theorem: 


ofot (r T ,x) = 2 


d 2 b T N(r T , b T ,x). 


(3.8) 


The photon wave function squared can be interpreted to be 

proportional to the probability for the 7* —>• qq splitting, where T and L refer 
to transverse and longitudinal polarizations of the virtual photon and / is the 
quark flavor. The transverse separation of the quarks is r 7 and the impact 
parameter in the collision is b The longitudinal momentum fraction of the 
photon carried by the quark is denoted by z. 

The photon to dipole splitting can be computed using quantum electro¬ 


dynamics written on the light cone as discussed in Sec. 2.2 The virtual 
photon wave functions squared (summed over spins and polarizations) can be 


obtained from the light cone calculations presented in Sec. |2.2| (see Eqs. (2.10) 


and (2.11)) and are [23J (see also my MSc thesis ra for a detailed derivation) 


o n 

[T*T]£ = a em e 2 f {[z 2 + (1 - z) 2 ]e 2 K 2 (er ) + m 2 f K$(er) 

7T 1 1 

1 / _ „2 n 2„2/! ; 


[T*T]£, = — -a em e 2 Q 2 z 2 (l - z) 2 K^(er), 
7r 1 


(3.9) 

(3.10) 


where e 2 = z( 1 — z)Q 2 + rnj and r = \vt |- 

Note that now all QCD dynamics is included in the virtual photon-proton 
cross section a 7 * p , and the “trivial” part of the process where the photon is 
emitted from the lepton is factorized out. The structure functions Fi and F 2 
can also be related directly to the virtual photon-proton cross sections, and 
in the high-energy limit we obtain m 


F2(x,Q 2 ) 
2 xFi(x, Q 2 ) 


_T_ a TP 

47T 2 a em tot 

Q 7 *p 

— 7 , - V T 

4vr 2 a em 


Q 2 

47T 2 « em 




(3.11) 

(3.12) 


Thus, F -2 measures the total cross section. At small-x the photon couples 
to the sea quarks that originate from gluons, and F 2 is related to the total 
gluonic density of the proton. The longitudinal structure function, on the 
other hand, measures the cross section for the longitudinally polarized virtual 
photon-proton scattering, as 


F l (x,Q 2 ) = F 2 (x,Q 2 ) - 2xF 1 (x,Q 2 ) 


Q 2 

47T Q?em 


7* 

L 


(3.13) 
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3.3 Extracting initial conditions for BK evolu¬ 
tion 


After the previous discussion we can compute the proton structure functions 
within the Color Glass Condensate framework given that we know the dipole 
amplitude N(rT,bx,x) for the dipole with size r t to scatter with impact 
parameter h?. The x (or energy) evolution of N is given by the Balitsky- 


Kovchegov equation (2.31) discussed in Sec. 2.4 However, the BK equation 


is a differential equation whose solution requires an initial condition. 

The initial condition for the BK equation, the dipole amplitude N(xt, &t, x) 
at some x = x$ is non-perturbative information. It can be obtained by 
parametrizing the amplitude and fitting the parameters to the experimental 
data. Note that this is exactly the same procedure that is used to determine 
the parton distribution functions: there one parametrizes the distributions at 
an initial scale Qq, and the distributions at higher scale Q 2 > Qq are obtained 
by solving the DGLAP equations. 

In current phenomenological literature the initial condition for the BK 


equation is often taken as an MV-model (see discussion in Sec. 2.3) inspired 
parametrization 


N(xt, x = xq) = 1 — exp 


(r 2 Q 


s,0 ) 


In 


r A, 


+ e c 


QCD 


(3.14) 


The fit parameters are Q 2 0 , the anomalous dimension 7 and an infrared cutoff 
e c . The parameter Q 2 0 can be roughly seen as a measure of the saturation 
scale, or as a gluon density at initial Bjorken-x. The anomalous dimension 7 
controls the behavior of the gluon distributions at large transverse momentum. 


Note that in Eq. (3.14) we did not include any impact parameter dependence. 


In practical calculations it is often assumed that one can replace 
2 / d 2 b T N(rT,b T ,x) -> a 0 N(r T ,x), 


(3.15) 


where cr 0 is twice the transverse area of the target, which must be fitted to 
the data. 

When the BK equation is solved with the running coupling corrections 


see Sec. 2.4.1), the strong coupling constant a s must be evaluated at some 


given distance scale. As discussed in Sec. 2.4.1 , the uncertainty in the Fourier 
transform of the expression of a s from the momentum space to the coordinate 
space is parametrized by factor C 2 , which is taken to be a fit parameter. See 


discussion in Sec. 2.4.2 for a discussion of a possible theoretical value for this 
parameter. 
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The first dipole amplitude fits with the running coupling BK equation to 
the HERA F 2 data, combined with CERN-SPS and FNAL measurements, 
were performed by the AAMS collaboration in Ref. m- This work was later 
improved by the AAMQS collaboration m where the combined measurements 
of the reduced cross section (see Eq. (3.5)) from the different HERA experi¬ 
ments 


(ZEUS and HI) with significantly smaller uncertainties were used. 
It was found in these fits that an anomalous dimension 7 > 1 and relatively 
large C 2 ~ 5... 10 (compared to theoretical calculation C 2 = e~ 27s ~ 0.3 
as discussed in Sec. 2.4.2[ ) were needed for a good description of the HERA 
data. The largish value of C 2 can be interpreted such that the leading order 
BK equation tends to give a too fast evolution speed compared to data, and 
this speed is reduced by a large C 2 . See also Sec. 2.5 for a discussion of NLO 
effects in the evolution speed of the BK equation. 

In paper m we studied a possibility to obtain a good fit to HERA data 
without introducing an anomalous dimension. A motivation for this was the 
fact that when an unintegrated (transverse momentum dependent) dipole 
gluon distribution Eq. (2.28) is computed from the dipole amplitude, one 
can obtain negative gluon densities using the initial condition obtained in 
Refs. [671 [79] before any BK evolution. Thus, we allowed the infrared cutoff 
in the modified MV model (3.14) to vary and parametrized it by introducing 
a fit parameter e c and fixing 7 = 1 . Another difference to the previous works 
was that we did not include any other experimental data except the combined 
HERA results, as thanks to the good precision of the HERA a r data, it 
anyway dominates the fit. 

We also only considered the light quarks instead of including charm as 
done in Ref. m, even though the charm contribution to the reduced cross 
section is measured by HERA |12j . The reason for this is that in Ref. |67] it 
was found that in order to obtain a good fit with heavy quarks one has to 
introduce additional parameters for the transverse area and initial saturation 
scale for the charm quarks. As we found this approach inconvenient, the 
heavy quarks were not included in the fit in paper m- 

We studied three different MV-model inspired parmetrizations for the 
initial condition. First, we used the standard MV model without modifications 
by fixing 7 = e c = 1. What is denoted by MV 7 is the same parametrization 
used also in previous works where the anomalous dimension 7 is a free 
parameter but we fix e c = 1. Similarly in the MV e parametrization we fixed 
7 = 1 but took the infrared cutoff parameter e c to be a fit parameter. The 
fit was performed to the HERA reduced cross section ay data in kinematical 
window x < 0.01 and Q 2 < 50 GeV 2 . I 11 this kinematics the CGC picture 
should be valid due to large gluon densities at small x, and the upper limit 
for Q 2 is needed as we did not include the Q 2 evolution via e.g. DGLAP 
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X 2 / d.o.f 

<32. o [GeV 2 ] 

Q'i [GeV 2 | 

7 

C 2 

e c 

? [mb] 

MV 

2.76 

0.104 

0.139 

1 

14.5 

1 

18.81 

MV 7 

1.17 

0.165 

0.245 

1.135 

6.35 

1 

16.45 

MV e 

1.15 

0.060 

0.238 

1 

7.2 

18.9 

16.36 


Table 3.1. Parameters from fits to HERA reduced cross section data for differ¬ 


ent initial conditions. Also the corresponding initial saturation scales Q 2 defined 
via Eq. (3.16) are shown. The parameters for the MV 7 initial condition are ob¬ 
tained by the AAMQS collaboration iU- 


equations. Note, however, that the results are not independent of Q 2 , as the 
virtuality of the photon sets the dominant transverse size r ~ 1/Q for the 
dipole. 

We observe that some modification 


The fit results are given in Table 3.1 


to the MV model is needed in order to get a good fit with y 2 /d.o.f ~ 1. For 
the MV 7 parametrization we obtain a very similar fit result as obtained by the 
AAMQS collaboration in Ref. [67], and small deviations can be explained by 
noticing the use of different data sets and small differences in the numerical 
setup. In our numerical calculations we will use the AAMQS result for the 


MV 7 parametrization, and we also show only their result in Table |3.1 

The initial saturation scale is controlled by the parameter Q S} o, but the 
numerical values for it are not directly comparable between the different 


parametrizations. Thus, we also show in Table |3T] the initial saturation scale, 
which is a physically relevant quantity, defined via the equation 


N (r 2 = ^ ) = 1 - e- 1 ' 2 . 


(3-16) 


The two parametrizations that give equally good fits to the data, MV 7 and 
MV e , also give a consistent result for the proton saturation scale at the initial 
x = 0.01: Q s « 0.5 GeV. 

The results obtained using the best fit values are compared with a part 


of the HERA reduced cross section data in Fig. 3.2 at different values of Q 2 
(note that Q 2 = 200 GeV" data is not included in the fit). Both modified 
parametrizations seem to give almost identical results that agree very well with 
the precise HERA data. The standard MV model is not too bad either, but 
the agreement is clearly worse than with the MV 7 or MV e parametrizations. 

In addition to the BK-evolved MV model parametrizations, similar fits have 
been done with other models for the dipole amplitudes. These are necessary 
especially if the impact parameter dependence of the dipole amplitude is 
needed, which is the case when dealing with e.g. diffractive scattering which 
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Figure 3.2. Reduced cross section o r computed using the different fitted 
parametrizations as an initial condition for the BK evolution compared with 
the combined HERA data from HI and ZEUS experiments [12] . Figure from 
paper (Vj. 


is the topic of Chapter |4j As an example we mention the IP sat model, 
which has an eikonlaized DGLAP-evolved gluon distribution function [80] fit 
in Refs. [78, 81] and the IIM model [82] that parametrizes the asymptotic 
properties of the BK equation, and has been fitted in Refs. |83j, S3]. A 
disadvantage of these parametrizations is that they do not include the full 
Bj or ken-a; evolution given by the BK equation. 


3.4 Dipole amplitude for a heavy nucleus 

Currently the DIS experiments with nuclei have only been done with fixed 
target experiments at SLAC [85] and by the NMC collaboration [861 1ST]. 
The kinematical range probed in these measurements is very limited for the 
small-x part of the phase space. This prevents us from performing a similar 
fit for the nuclei as what is done for the protons. 

In the coming decades there are plans to build next-generation colliders 
that could study lepton-nucleus scattering at high energy. In the United 
States, the plans are to build an Electron-Ion Collider (EIC) [88] by adding 
an electron beam to the BNL-RHIC accelerator, which currently collides 
protons and different nuclei, or adding a nuclear beam to the JLAB-CEBAF 
facility. At CERN, the possibility to add an electron beam to the Large 
Hadron Collider, which can accelerate protons and lead nuclei, is considered 
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as a way to build an LHeC collider [89j . 

Before experimental small-rr nuclear DIS data becomes available another 
approach must be taken in order to obtain the dipole-nucleus scattering 
amplitude N A . We use the optical Glauber model (for a review, see e.g. 
Ref. [90]) to generalize dipole-proton scattering amplitude to dipole-nucleus 
scattering to obtain an initial condition for the BK evolution. We presented 
this method in paper [V] . 

Let us first write the total dipole-proton cross section for the dipole with 
a transverse separation r as 


ct£j = CTo N(r), 


(3.17) 


where <7o is the result of the impact parameter integral in Eq. (3.8), and was 


fitted to HERA data in Sec. 3.3 In the dilute limit where the dipole is very 


small, the dipole-nucleus cross section should be just A times dipole-proton 
cross section (where A is the mass number of the nucleus). On the other 
hand, for large dipoles we should have 


^ 2N(r T M) -2. 


(3.18) 


These requirements are satisfied with an exponentiated dipole-nucleus scat¬ 
tering amplitude 


N a (vt, hr) 


1 — exp 


AT A (b T ) 



(3.19) 


where T A is the transverse density of the nucleus at impact parameter br 
which, in practice, is obtained by integrating the nuclear Woods-Saxon EB 
density distribution over the longitudinal direction. This form corresponds 
to an average of the dipole-nucleus scattering amplitude over the fluctuating 
positions of the nucleons in the nucleus, see e.g. Refs. [ 80] 9_2 . 

The scattering matrix S A = 1 — N A calculated from Eq. 
a limiting value 

exp ~ ^ 1/3) ) (3.20) 


(3.19) approaches 


instead of zero, in the large dipole limit. This will cause unphysical oscillations 
to the gluon distribution computed from the dipole amplitude (see discussion 


in Sec. 2.3). Therefore we expand the dipole-proton cross section crL as 


T P - 

>qq 


TV 


00 - 






In 


|^t|Aqci 


+ e c • e 


(3.21) 
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N(r 2 = 2/Ql) = 1 - e - 1 / 2 



Figure 3.3. Saturation scale for the 
proton and the lead nucleus at central 
and semi-central collisions as a func¬ 
tion of x calculated using the MV 7 
parametrization. 


N(r 2 = 2/Q 2 ) = 1 - e- 1 / 2 



Figure 3.4. Impact parameter depen¬ 
dence of the saturation scale of the lead 
nucleus at two different values of x. The 
dashed lines correspond to the proton 
saturation scales. Figure from paper 0 - 


Combining this with Eq. (3.19) we obtain the dipole-nuclues amplitude at 
x = X 0 : 


N A (r T , b T ) = 1 — exp 


-AT A {b T ) 


(To {TtQI ,0 


In 


|Tt|Aqc 


+ e c 


D 


(3.22) 

which is our result, from paper (V], for the dipole-nucleus amplitude N A . One 
advantage of the MV e parametrization, where 7 = 1, can now been seen from 
Eq. (3.22): in this case there is no ambiguity whether the initial saturation 
scale should scale as Ta or . 

To obtain the energy (or Bjorken-x) evolution of the dipole-nucleus ampli¬ 
tude we use the obtained fit parameters from lepton-proton DIS quoted in 
Table 3T and solve the BK equation separately at different impact parameters 
using Eq. (3.22) as an initial condition. Note that in principle one should 
evolve the initial condition using an impact parameter dependent BK equation 


which is not currently possible in practice, as discussed in Sec. 2.4.1 


The obtained saturation scales for the proton and the lead nucleus are 
shown in Figs. |3.3 and T4 For the lead the saturation scale is shown for 
the central and mid-central impact parameters (corresponding to 0 — 20% 
and 40 — 60% most central collisions, defined from the Optical Glauber 
picture |90|). The saturation scale of the nucleus is significantly larger than 
that of the proton still in mid-central centrality classes, and falls below the 
proton saturation scale only at b > 6 fm, which corresponds to centrality 
classes > 70%. In that region our parametrization is not exactly valid any 
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more, as the BIv evolution would make the nucleus to grow rapidly from 
the dilute edges. In phenomenological applications we choose to use the 
parametrization (3.22) only in the region where the saturation scale in the 
nucleus is larger than in the proton. 
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Chapter 4 

Exclusive vector meson 
production 


4.1 Diffraction in scattering experiments 


In high-energy scattering experiments diffractive events are experimentally 
defined such that a large rapidity gap is present in the event. A rapidity gap 
means that there is an interval in rapidity (a few units wide) where there are 
no particles produced. Elastic scattering is a simple example of a diffractive 
event, but it is also possible to have other events where, for example, a single 
particle is produced at central rapidity and the scattered particles continue 
at very forward/backward rapidities. 

An example of a diffractive scattering process is shown schematically in 


Fig. 4.1 where in a lepton-proton scattering the final state particles are the 


scattered lepton scattered proton P' and the produced system X. We 
focus on exclusive vector meson production, where the produced system X is 
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for example a J/f meson. The scattered proton can either remain intact or 
break np, but there can not be exchange of color charge between the proton 
and the scattered system, or otherwise the breakup of the color strings would 
fill the rapidity space between the proton and the system X. In addition 
to lepton-proton scattering, diffractive scattering can be studied also with 
nuclear targets and in hadronic collisions (see Sec. 4.5). 

The photon-proton interaction can be described in terms of an exchange 


of a color neutral object, known as the pomeron. It is represented in Fig. 4.1 


by a zigzag line. The kinematical variable xp has the similar interpretation as 
Bjorken-x in DIS: it describes the fraction of the incoming proton (or nucleon) 
longitudinal momentum carried by the pomeron: 

(P — P') ■ q M 2 + Q 2 — t 


x P = 


P ■ q 


2 ’ 


W 2 + Q 2 - rn N 


(4.1) 


where q is the four-momentum of the emitted photon. In addition, in order 
to describe the kinematics of diffractive deep inelastic scattering (DDIS) 
processes, we define the momentum transfer t — (P — P') 2 and use the 
virtuality of the photon Q 2 = — q 2 and the center-of-mass energy of the 
photon-nucleon system W 2 = (P + q) 2 . 

Diffractive scattering processes can be used to probe the target gluon 
distribution very efficiently. This is because at leading order there is an 
exchange of two gluons between the dipole and the target in the scattering 
amplitude, as there can not be net exchange of color charge. Thus, the 
amplitude is proportional to the leading order gluon distribution function 
xg, and the total cross section behaves like gluon distribution squared, see 
Ref. [93]. In addition to gluon densities, exclusive vector meson production 
can be used to probe the spatial distribution of gluons within the proton or a 
nucleus, see e.g. Refs. [941496] . 

Currently the most important diffractive DIS measurements come from 
HERA, where diffraction off a proton was studied. For nuclear targets, 
diffractive DIS has so far been measured only in fixed target experiments 
such as E665 ra and NMC [H8] at relatively low center-of-mass energies. In 
the future, if an electron-ion collider (EIC) [SB] or LHeC [S9J is realized, it 
would open a new era for studies of diffractive DIS off nuclei. 


4.2 Exclusive vector meson production in the 
dipole picture 

Let us consider diffractive vector meson production in high-energy deep 
inelastic scattering, where the target can be either a proton or a nucleus. 
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Figure 4.2. Diffractive scattering in dipole picture. 


The theoretical framework for description of diffractive scattering events, 
the so called Good-Walker picture, was developed in Ref. (99J. In this 
framework one has to find the states that diagonalize the imaginary part 
of the scattering T-matrix, which at high energy are the ones where the 
virtual photon fluctuates into a quark-antiquark dipole long before the dipole 
interacts with the target. The dipole interacts elastically with the target and 
finally forms a vector meson (p, u, J /T, T,...). This is shown schematically 
in Fig. |4.2| Note that as we consider here high-energy processes where the 


eikonal approximation is valid, the transverse positions of the quarks can be 
taken to be fixed during the interaction. 

In order to calculate the cross section one has to compute virtual pho¬ 
ton splitting to a quark-antiquark dipole (the virtual photon wave function 
which was already discussed in Sec. 3.2. In addition, the formation 
of the vector meson requires knowledge of the vector meson wave function 
T\/, which is in general non-perturbative. In practice for heavy mesons a 
large contribution to the cross section will come from distance scales set by 
the meson mass \rx\ ~ 1/m. The elastic dipole-target scattering is described 
by the dipole amplitude N. 

The cross section for the diffractive vector meson production is derived in 
Refs. [93, HUD ] (see also Refs. (2UIZ8I). The scattering amplitude is 

U ' 4 - 1 ’ - J dV J I'M Or. z) J 

(4.2) 

and the differential cross section can be computed as 


da 7 


*A->VA 


di 


167T 


\A 


7*A->VA| 2 


(4.3) 


Here z is the longitudinal momentum fraction of the virtual photon carried 
by the quark, bx is the impact parameter and A = \f—t is the transverse 
momentum transfer. 

















The scattering amplitude (4.2) has a clear physical interpretation: first the 
virtual photon splits into quark-antiquark dipole with probability amplitude 
and the dipole then scatters elastically off the target color field 
described by the dipole amplitude N. Finally, the quark-antiquark pair forms 
the vector meson with probability amplitude T. The momentum transfer 
dependence of the scattering amplitude originates from the Fourier transform 
of the dipole amplitude. It is then necessary to have a detailed description 
of the impact parameter profile of the target, in contrast to the DIS cross 
section calculations where we only had to calculate the integral over the 
impact parameter space in Chapter [3j 

As the vector meson wave function can not be computed from first princi¬ 
ples using perturbative techniques, different models for it exist in the literature. 
The wave function is obtained by assuming that the vector meson is pre¬ 
dominantly a quark-antiquark state that has the same spin and polarization 
structure as the virtual photon. The model parameters are fixed by requiring 
the wave function to reproduce the measured decay width to the electron 
channel (which turns out to be proportional to the wave function at origin) 
and that it is correctly normalized. In this work we study the dependence of 
our results on the model uncertainties by using two different models for the 
vector meson wave function, called “Gaus-LC” and “boosted Gaussian”. These 
wave functions are described in more detail and fitted to data in Ref. |78j . 

In order to calculate the evolution of the dipole amplitude with impact 
parameter dependence a solution to the impact parameter dependent BK 
equation would be needed. However, as discussed in Sec. |3.4[ it is not 
currently available. Instead of a BK evolved dipole amplitude we will use 
different parametrizations for N(xT,bT,Xp) which include impact parameter 
dependence. 

We use mostly the so called IPsat model [HD], where the dipole amplitude 
is obtained from an eikonalized DGLAP-evolved gluon distribution function 
xg whose initial condition is fitted to the HERA deep inelastic scattering data 
in Ref. The transverse profile of the proton is assumed to be Gaussian, 
and the width of the distribution is a fit parameter. 

The second dipole amplitude that we use here to study the dependence 
of our results on the details of the dipole amplitude is the so called IIM 
model [82], which is a parametrization including the most important features 
o the BK evolution. The model parameters are again fitted to the HERA 
data-in Ref. [83] . The calculated diffractive J/'F production cross section in 


1 After the completion of paper (Ij, where the parametrizations were used, newer fits for 
both IIM and IPsat models to the much more precise combined HERA data have been 

published mm 
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electron-proton scattering is shown in Fig. 4.3 as a function of Q 2 , where the 
results are compared with the HI [ 101 ] and ZEUS [ 102J data. The “Factorized 
IPsat” refers to a parametrization where the impact parameter dependence is 


factorized from the IPsat model, see discussion in Sec. 4.3 


4.3 Diffraction off a nucleus 


When diffractive scattering with nuclear targets is considered, the events can 
be divided into two classes. In coherent diffraction the target nucleus remains 
fully intact, whereas in incoherent diffraction the p t kick given to the nucleus 
is large enough to break it up still preserving the rapidity gap. 

As shown in Ref. [93], the coherent cross section corresponds to performing 
the average over the nucleon configurations (or the nuclear wave function) at 
the scattering amplitude level, and squaring the amplitude gives the cross 
section. Similarly, averaging the cross section instead of the amplitude over 
the nucleon positions gives the sum of incoherent and coherent cross section, 
called quasielastic cross section. 

Coherent diffraction dominates at small momentum transfer t ~ — 1 /R\ 
(where R A is the nuclear radius), where the dipole scatters coherently off 
the whole nucleus. At larger momentum transfer t ~ — 1 /R v (where R p is 
the proton radius) the probed objects are individual nucleons instead of the 
nucleus. Experimentally, the t dependence of the coherent cross section is a 
challenging measurement due to the difficulties of measuring small momentum 
transfer and the intactness of the nucleus. 

The cross section for the quasielastic vector meson production can be 
written as (see e.g. Ref. [75] ) 


da~'* A ^vA 

d t 


R 2 9 ( i + l 2 ) 

167T 


(|3.(xp, Q 2 , At| 2 )at, 


where the average over the nucleon configurations is denoted by 


(4.4) 


<0({frr,i})> 



[d 2 b Tii T A {h Tti )\ 0({b T>i }). 


(4.5) 


The nucleon positions bx,i are assumed to be independent, and as a nuclear 
thickness function T A we use the Woods-Saxon distribution f91j. The coherent 
cross section is obtained by averaging the amplitude before squaring it, | (A)n | 2 , 
and the incoherent cross section is variance (|*4.| 2 )tv — I(*4)jv| 2 that measures the 
fluctuations of the gluon density inside the nucleus. The averaged amplitude 
(3 )tv is a smooth function of the impact parameter br, and its Fourier transfer 
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vanishes rapidly at A > I / R A . Therefore, at large momentum transfer A, 
the quasielastic cross section ( |4.4 ) is almost purely incoherent. 

The factor 1 + (3 2 takes into account the correction from the real part of 
the scattering amplitude and R g corrects for the so called skewedness effect, 
i.e. that the gluons in the target are probed at different x, see Ref. [103j. The 
corrections are calculated following the prescription of Ref. [ 104 ]. In order 
to derive the cross section for the incoherent cross section we modify the 
IPsat model and use a factorized approximation for the scattering matrix 
S P {r T , b T , xp) = l — N p (r T , b T , x P ) as 


S p (r T , b T , zp) = 1 — T p (b T )N p (r T , x P ). (4.6) 

In paper [Ij it was shown that this approximation changes the cross section 
very little when calculating diffractive vector meson production with proton 
targets. 

The cross section for the coherent diffractive vector meson production is 
derived in Ref. [80], the result being 


(A(x ¥ ,Q 2 , A t ))n = f ^d 2 r T d 2 6 T e-^ A -[^t,^^](r T ,g 2 ,z) 

x 2 [1 - exp {—27 t B p AT A (b T )N(r T , z P )}]. (4.7) 

Here B p is the width of the Gaussian density distribution of the proton. The 
quasielastic (which is almost purely incoherent at large \t\) cross section was 
derived in paper [I] to be 


(|^4(zp, Q 2 , At) | )n — 
x [T^ 7 


167 rB p A [ d 2 b T d 2 r T di 2 r' T ^f-^f- 
J Air 47 t 

Q 2 , ^)[Tt,T 7 *-^](4, Q 2 , zf) 


x g-BpA?, e -2nB p AT A (b T )[N(r T ,x v )+N(r' T ,x P )} 


7T BpN(r T , x v )N(r' r , x v )T A {b T ) 

1 - 2nB p T A (b T ) [N(r T , z P ) + N(r ' T , z P )]' 


(4.8) 


The squared amplitude is proportional to A times the squared dipole-proton 
amplitude, corresponding to independent scattering off the nucleons. This is 
multiplied by a nuclear attenuation factor 


e -2TrB p AT A (b T )[N(r T ,x r )+N(r' T ,x r )\ 

1 - 27 vB p T A (b) [N(r T , x P ) + N(r' T , z P )] 


^ e -^{A-l)T A {b T )[N{r T ,xt)+N{r' T ,x r )\ 

(4.9) 
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7 *p—tJ/'ffp, W =90 GeV 



Figure 4.3. Diffractive J/'I' produc¬ 
tion cross section as a function of Q 2 
compared to different dipole model cal¬ 
culations. Figure from paper [I]. 


A= 197, Q 2 =0 GeV 2 , i p =0.001 



Figure 4.4. The quasielastic and coher¬ 
ent diffractive J/’k cross sections in gold 
nuclei. Figure from paper [IJ. 


which accounts for the requirement that the dipole can not scatter inelastically 
off the remaining A — 1 nucleons in order to keep the event diffractive, see 
paper p] for details. 

The calculated incoherent and coherent diffractive J/T photo prod uction 
( Q 2 = 0 GeV 2 ) cross sections in 7 *-gold scattering are shown in Fig. 4.4 where 


the differential cross section is shown as a function of squared transverse 
momentum transfer t. Note that our approximation for the incoherent cross 
section is not valid at small \t\, and the corresponding part of the distribution 
has been left out. To estimate the saturation effects we also calculate the 
cross section using the “IPnonsat” model, where the IPsat model is linearized 
such that the dipole amplitude is proportional to at all dipole sizes. When 
this parametrization is used, the incoherent cross section is explicitly A times 
the dipole-proton cross section (as calculated e.g. in Ref. pH]). The two 
saturation model parametrizations (IPsat and IIM) give comparable results 
for both coherent and incoherent cross sections, and the incoherent cross 
section is strongly suppressed by the saturation effects. 

The origin of the suppression can be seen from the nuclear attenuation 
factor (4.9): in the black disk limit N(rx,xp) = 1 (which is the case when 
r T > 1 /Qs), and the nuclear attenuation factor at small impact parameters 
behaves like ~ e -° xa 1 / 3 ^ w jq c p can b e seen by noticing that Ta^t) ~ A '~ 2 / 3 
when Ifonl is not very large. Thus the contribution from the center of the 
nucleus to the incoherent cross section is suppressed, and only the scattering 
from the edges of the nucleus contributes. As the edge of the nucleus has an 
area ~ 2 TtRAd (where d is the thickness of the edge), the incoherent cross 
section behaves as ~ A 1 / 3 in the black disk limit. On the other hand, in the 
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A = 197, |*| =0.5 GeV 2 



Figure 4.5. Nuclear suppression factor 
of the cross sections as a function of 
Q 2 for diffractive J/’F production. The 
upper thin curves correspond to xp = 
1CT 2 , and the lower thick curves are 
calculated at xp = 1CF 4 . Figure from 
paper [I]. 


A = 197, |i| =0.5 GeV 2 



Figure 4.6. Nuclear suppression factor 
of the cross sections as a function of 
xp for diffractive J/\F production at 
photoproduction and at Q 2 = 10 GeV 2 . 
Figure from paper [I], 


dilute limit (with no saturation effects) the attenuation factor goes to unity 
as the dipole amplitude is small, and the cross section is proportional to A. 

To quantify the nuclear effects we calculate the nuclear suppression factor 
defined as the ration of the 7 *-nucleus and 7 *-proton cross sections normalized 
by the number of nucleons A in the gold nucleus. The scale Q 2 dependence 
of the suppression factor is shown in Fig. 4.5 and the same quantity as a 
function of x> is plotted in Fig. 4.6 The suppression is larger at small Q 2 
and at smaller xp where the saturation effects should be large. Note that 
the dominant dipole size is r ~ 1/Q, and thus at large Q 2 a dilute region 
of the nucleus is probed where no nuclear effects are expected. Similarly 
at smaller xp the non-linear effects should be larger, which can be clearly 
seen from Fig. |4.6| where the suppression grows as xp decreases. The model 
uncertainties are quantified by calculating the ratio with both IPsat and IIM 
parametrizations, the difference being largest at small Q 2 and xp. 

In Fig. 4J3 the results are calculated also using a nonsaturated dipole- 
nucleon cross section (“IPsat, nonsatp”), which corresponds to including 
unitarity (or equivalently saturation) effects at the level of the nucleus but not 
for a single nucleon. With this parametrization a much larger suppression is 
obtained than with the IPsat model. This shows that the nuclear suppression 
ratio is also sensitive to the saturation effects at the proton level. 

The fact that the different parametrizations (IPsat and IIM) differ signifi¬ 
cantly when calculating spectra or nuclear suppression factor shows that the 
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diffractive electron-nucleus scattering has a discriminating power to separate 
between different dipole models, and it can be used to constrain the model 
uncertainties. Note that here both parametrizations are fitted to HERA DIS 
dataQ but the differences in eA scattering can be up to 50%. The calculations 
presented here are done in the kinematical region that would be available in 
a future electron-ion collider. 


4.4 Centrality in diffractive events at an Elec¬ 
tron Ion Collider 

Centrality selection based on event multiplicity has been a powerful tool 
to learn about the QCD dynamics in heavy ion collisions. The multiplicity 
selection has been applied also in proton-nucleus and proton-proton colli¬ 
sions. Recent studies of high multiplicity proton-proton and proton-nucleus 
events have revealed many interesting phenomena such as long range angular 
correlations (see e.g. Ref. [9]) that are not always visible in minimum bias 
events. 

In the previous section the diffractive vector meson production cross 
section was calculated for lepton (or photon)-nucleus scattering. It was also 
found that coherent diffraction dominates at small \t\, whereas incoherent 
diffraction is the only component at \t\ ~ 1 /R p where the p? kick is localized 
on an area comparable to the nucleon area. Incoherent diffraction then probes 
fluctuations at the distance scale of the order of the nucleon size pn HUH], 
and the incoherent cross section can be expected to depend on the impact 
parameter of the photon-nucleus collision^} Triggering on the most central 
collisions we could probe the densest region of the nucleus where the saturation 
scale Q s is expected to be enhanced relative to the minimum bias events. 

In heavy ion collisions the centrality selection can be done relatively easily 
according to the total multiplicity of the event, such that the events with 
largest multiplicity are the ones where the two nuclei collided with zero 
impact parameter. The multiplicity classes in proton-nucleus collisions are 
discussed later in Chapter [5j What we propose in paper ivq is to use the 
proton multiplicities in the “Roman pot” detectors as a measure of the impact 
parameter in incoherent diffraction. 

Consider a nucleus participating in a event where a vector meson is 

2 Although we have made a simplifying assumption for the impact parameter dependence 
of the IPsat model. 

3 With \t\ > 1 /R p the cross section would also be sensitive to sub-nucleon scale fluctua¬ 
tions which are not completely included in this work. 
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7 + A -> V 1 /V 2 + A*, x r = 0.005 



Figure 4.7. EIC acceptance for pro¬ 
tons, as a function of pr, in a Roman 
pot detector. Figure from Ref. [Mi- 



Figure 4.8. Ratio of two incoherent 
vector meson production cross sections 
in central events relative to minimum 
bias as a function of Q 2 . Figure from 
paper ED 


produced diffractively. If a nucleon in the nucleus receives a large momentum 
kick ~ \f\t\, it escapes from the nucleus and can scatter off other nucleons 
on its path out. When the scattering takes place close to the center of the 
nucleus, there are more nucleons off which the participating nucleon can 
scatter and on average we expect more “ballistic” nucleons to be produced 
than in peripheral events. These nucleons travel along the beam pipe but 
their trajectories differ slightly from that of the original nucleus and can 
be measured in the Roman pot detectors located in the beam pipe outside 
the main detectors. The simulated acceptance of the EIC [SB] Roman pot 


detector is shown in Fig. 4.7 where it can be seen that the acceptance is good 
(generated and accepted spectra are the same) in the region of pr that the 
ballistic protons would have [^] Note that now we can only discuss about 
ballistic protons, as the neutrons are not bent by the magnetic held and do 
not reach the Roman pots. 

In diffractive events the nucleus is usually left in an excited state. When 
the nucleus returns to its ground state nucleons can evaporate according to 
a thermal spectrum. In the laboratory frame these nucleons should have 
significantly smaller transverse momentum than 400 MeV... 1 GeV that we 
expect for many ballistic protons, and are not detected in the Roman pots. 

The produced neutrons from both of these channels can be measured 
in Zero Degree Calorimeters (ZDC), but there the separation between the 


4 The simulation is done for Deeply Virtual Compton Scattering in e+p collision, in e+A 
scattering the acceptance will shift to lower pr due to the higher magnetic fields required 
by the different mass-to-charge ratio. 
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ballistic and thermal components is difficult. Thus we do not expect the 
ZDC measurements to give sufficient information to divide diffractive events 
into the centrality classes. Note, however, that in other processes such as 
single inclusive multiplicities and dihadron correlations the ZDC energy has 
a potential to be a measure of the event centrality mm. 

As a model example we study the centrality dependence of the incoher¬ 
ent diffractive vector meson production cross section using the framework 
described earlier in this Chapter. “Central” events, where the proton multi¬ 
plicities in the Roman pot should on average be large, are defined by setting 
the impact parameter to zero when calculating the incoherent cross section 
from Eq. (4.8). The production cross section is compared with the minimum 


bias cross sections obtained by integrating over all impact parameters. To 
cancel the theoretical uncertainties associated with e.g. the vector meson 
wave functions, we consider the Q 2 dependence of the double ratio 


a( 7 *A V\A*) ja(YA V 2 A*) 


central 


cr(y*A —> V\A*) /< 7 ( 7 *A —>• V 2 A*) 


(4.10) 


minimum bias 


where V\ and V 2 refer to different vector meson species. In this analysis, we 
include J/T, p and 0. In our approximation the real part and skewedness 
corrections and t dependence cancel in the double ratio. 


The results are shown in Fig. |4.8| where the double ratio is calculated at 
xp = 0.005 which is within the EIC kinematics. A significant enhancement in 
the J/T/p and J/T /0 ratios is seen at small Q 2 . This enhancement can be 
understood as follows. As J/T is more massive, its wave function is peaked 
at smaller dipoles and the cross section is dominated by configurations with 
r 2 Q 2 <C 1 also at low Q 2 . On the other hand for larger p and 0 mesons the 
typical configurations have r 2 Q 2 > 1 at the same Q s , and the saturation 
effects affect the cross section more in central events (where the saturation 
scale Q s is larger) than in the minimum bias events when Q 2 is decreased, 
which increases the double ratio. 

Another way to think of this is to note that as r 2 Q 2 > 1 in the production 
of larger mesons (p and 0), the cross sections in central and minimum bias 
events are proportional to the geometric area of the interaction and cancel 
from the double ratio. Moreover, as the J/T production cross section goes like 
Q 4 (as the dipole amplitude is proportional to Q 2 in the dilute region), the 
double ratio becomes just the ratio of the saturation scales Q 4 rpr , tra i /Qj 1Tlir , hias . 
At large enough Q 2 also p and 0 cross sections become dominated by dipole 
sizes where r 2 Q 2 <C 1 and the numerators and denominators cancel separately 
from the double ratio, and we obtain unity as can be seen from Fig. |4.8| Since 
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p and 0 have roughly the same size, they are simultaneously probing either 
dilute or dense limits ( r 2 Q 2 s <C 1 or r 2 Q 2 > 1) at fixed Q 2 , and thus almost 
no Q 2 evolution can be seen in the 0/p double ratio. 

The calculation shown in Fig. |4.8| is a direct measure of the nuclear 
enhancement originating from the non-linear gluon dynamics described in 
the CGC framework. Even though we can not tell how many protons in 
the Roman pot corresponds to very central events, experimentally the most 
central event class can be defined with events that have largest number of 
ballistic protons in the Roman pot. These events, on average, are more 
central than the minimum bias results. We expect that the results for the 
double ratio will be very different when calculating in models where the 
non-linear QCD dynamics is included differently than in the dipole model. 
Thus, even though quantitative data comparisons require more work (e.g. 
combining the presented calculation with a more detailed model for the 
nuclear breakup), interesting qualitative results can be obtained from the first 
EIC measurements. 


4.5 Ultraperipheral heavy ion collisions 

As discussed in this Chapter, diffraction in lepton-nucleus collision has a 
potential to be a powerful tool to study the properties of the dense QCD 
matter, that is, the small-x structure of nuclei and protons. In deep inelastic 
scattering and lepton-nucleus diffraction the QCD dynamics is encoded in 
the virtual photon-nucleus scattering and the role of the lepton is to act as a 
source of virtual photons. 

An electron-nucleus collider would be an ideal tool to study, for example, 
exclusive vector meson production, but new experimental facilities such as the 
EIC |88] or the LHeC |89j colliders are needed for these studies. With currently 
available experiments, one possibility to study photon-nucleus scattering is 
to consider scattering processes where two hadrons collide with such a large 
impact parameter that the strong interactions are suppressed and can be 
neglected. Then, the role of one of the colliding (charged) objects is to act 
as a moving electromagnetic charge emitting virtual photons that scatter off 
the second hadron. This is exactly what is done in ultraperipheral heavy ion 
collisions. 

Let us consider a scattering of two heavy nuclei at high energy such that 
the impact parameter is larger than twice the nuclear radius, b > 2 R A , as 
shown schematically in Fig. |4.9| In these events it is possible to have photon- 
nucleus or photon-photon scattering. The later is suppressed by additional 
powers of electromagnetic coupling a em and is not considered here. 
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Figure 4.9. Two heavy nuclei colliding with large impact parameter b at high 
energy. 


Following the review [ 108 J the cross section for ultraperipheral hadron- 
hadron collision a AA can be written as 


a = dui 


n(oj) 


-a 


7 A 


UJ 


(4.11) 


where a lA is the photon-hadron cross section. The photon flux integrated 
over all possible impact parameters b > b m \ n = 2 Ra is 


n{uj) 


2Z 2 oi em 

7 t/3 




K) - A'oK)) 


(4.12) 


where £, = 2uRa/('JI3), Z is the number of protons in the nucleus and 7 is 
the Lorentz boost factor of the beam in the center of mass frame and ~ 1 
is the velocity of the incoming hadron. 

Let us consider the production of a vector meson with mass My, momen¬ 
tum pv and rapidity y in the laboratory frame in a process where a photon 
(momentum q, energy cu) is emitted from the nucleus that is moving along 
the z axis in the positive direction, and the nucleon it is scattering off has a 
momentum P. Assuming that the transverse momentum of the diffractively 
produced meson is small compared to its mass we obtain 


Py = —=m T e y ~ —=M v e y , 


V2 


V2 


(4.13) 
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where nir = \JM 2 + pj, is the transverse mass. On the other hand, momen¬ 
tum conservation gives 


p+ = q + + X P + = 


(4.14) 


as P + ~ 0, q + = y/2 uj and x is the longitudinal momentum transfer from the 
second nucleus. This gives 


U! = 



(4.15) 


Differentiating Eq. (4.11) with respect to to we find 


da AA _ n(u) A 
dev ~ u * 

and noticing that duj/uj = dy we can write 


da 

dy 


J d tn(uj) 


dt 


(4.16) 


(4.17) 


This allows us to calculate the differential cross sections as we know how to 
compute da lA /dt, as discussed in previous Sections. Note that experimentally 
one can not determine which one of the nuclei acted as a source of photons 
and from which nucleus the photon scattered off. Thus, both processes must 
be taken into account by calculating 


dcr 

dy 


n{y)o lAl (y) +n(-y)ti lA, (-y). 


(4.18) 


Let us then limit ourselves into the photoproduction region where Q 2 = 0 
and calculate the kinematical invariant xp. First, we observe that by definition 
W 2 = (q + P ) 2 « 2q ■ P « 2 q + P~ = 2a j^/s, as q + = \/2u and P~ = -^=1E 
with y/s = 2 E where E is the energy of the nucleon. Now using the definition 


of aip from Eq. (4.1) gives 


x P 


M 2 


Mi 


V -y 

e y . 


(4.19) 


Because scattering off both nuclei is possible, in vector meson production at 
large rapidities both the small and large x structure of the nucleus is probed: 
a small -x photon can scatter off a large-x gluon or vice versa. The dipole 
model calculations are only valid at small xp, which limits the applicability 
of the model close to central rapidity. 

The ALICE collaboration has measured both coherent and incoherent 
diffractive J/4/ production |109| differentially in rapidity. Thanks to the large 
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Figure 4.10. Coherent diffractive J/T 
production cross section computed us¬ 
ing the factorized IPsat (“flPsat”) and 
IIM parametrizations compared with 
the ALICE data [ 109) . Thin lines are 
obtained by using the Boosted Gaussian 
wave function for the J/’F, and for the 
thick hens Gaus-LC is used. Figure from 
paper El- 


Figure 4.11. Incoherent diffractive 
J/T production calculated in dipole 
model (“LM-fIPsat”) and using a leading 
twist approximation (“RSC-LTA”) com¬ 
pared to ALICE data [109] and to the 
STARLIGHT event generator. Figure 
from Ref. m- 


center-of-mass energy, the dipole model results are valid within \y\ < 2 for 
the J/\F production, as in this region xp < 0.02. In case of RHIC the smaller 
energy used in the experiment limits the applicability of the framework such 
that the production cross section at central rapidity can barely be computed, 
but no rapidity dependence is obtainable. 

In Papers [T] and BE! we compute the diffractive J/\\> production cross 
section using the same dipole models (IPsat and IIM parametrizations) as 
discussed earlier in this Chapter. Recall that the input to these calculations 
comes from electron-proton DIS measurements combined with standard nu¬ 
clear geometry. First, we compare our results for coherent diffractive J/'F 
production in ultraperipheral lead-lead collisions at yj snn = 2.76 TeV to 
the ALICE data [109] . The results are shown in Fig. 4.10, where the model 


uncertainty is quantified by calculating the cross section using both IPsat 
and IIM parametrizations and using two different wave functions for the 7 * 
and J/ X I> overall (Boosted Gaussian and Gaus-LC from Ref. [78]). 

Our results slightly overshoot the data, but the rapidity dependence comes 
out correctly, which means that evolution of the dipole cross section as a 
function of x-p is roughly correct. All parametrizations give consistently 
&a/(\y\y = o /da/dy y= 2 = 1.41 ... 1.46, thus the prediction for the rapidity 
dependence is more robust. The absolute normalization is also more strongly 
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model dependent that the shape of the distribution. It is also important to 
note that the two different wave functions used here have largest difference at 

Q 2 = 0. 


There are many calculations of the coherent vector meson production cross 
sections in the market (for references and comparison to the ALICE data see 
Ref. [1091]). but the situation is completely different for incoherent diffraction. 
In Fig. 4.11 we show the ALICE result for the incoherent diffractive J/T 
photoproduction cross section in ultraperipheral collisions. The dipole model 
predictions from paper IE! shown in the figure (“LM-fIPsat”) is obtained by 
using the IPsat model and Gaus-LC wave function. If the Boosted Gaussian 
wave function is used, the normalization changes and a good agreement with 
the ALICE data is obtained, but simultaneously the coherent cross section is 
more overestimated. 

PHENIX collaboration at RHIC has also measured coherent diffractive 
J/tf production cross section at midrapidity obtaining 76±34/rb |110j . Using 
the same parametrizations as when calculating the incoherent cross section in 
Fig. 4.11 we obtain 109/ib, being consistent with the data given the relatively 
large experimental uncertainty. 

As shown above, the absolute normalization depends strongly on the 
details of the model parametrizations, but the rapidity dependence is a more 
solid CGC prediction. Simultaneous comparison of the model calculations to 
experimentally measured coherent and incoherent vector meson production 
cross sections has a clear potential to constrain the dipole model uncertain¬ 
ties. In addition to ,//T, it is possible to also consider production of other 
mesons such as H/(2S') that are experimentally measured by ALICE [lTl| and 
calculated using our framework m- 
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Chapter 5 


Single inclusive particle 
production 

5.1 Particle spectra 

One of the simplest collider experiments are measurements of the single 
particle spectra. The total multiplicities in proton-proton and proton-nucleus 
collisions are dominated by low-p^ ~ Aqcd particles that originate from 
partons with low transverse momentum pr whose formation can not be 
described using perturbative techniques. When the particle production is 
measured differentially in we can limit ourselves into the kinematical 
region where perturbative techniques can be applied. 

In collider experiments hadrons, not individual quarks or gluons, are 
collided. Perturbative calculations, on the other hand, can be used to describe 
only the partonic scattering processes. The partonic density of the proton is 
given in terms of the parton distribution function (PDF) xfi(x, Q 2 ), which 
gives the probability to pick a parton i from the proton with momentum 
fraction x at scale Q 2 . The scale evolution of / is given in terms of the DGLAP 
equations, and the initial condition is obtained from fits to the HERA deep 
inelastic scattering data. In this thesis the CTEQ parton distribution function 
set cni is used. 

The partons are not measured in detectors, as the produced parton 
hadronizes to a cascade of color neutral particles. This transition is described 
by a fragmentation function D^h(z,Q 2 ), which gives the probability to 
produce a hadron h which carries the fraction z of the parent parton k 
momentum, at scale Q 2 . The scale evolution of IJk^h is again given by the 
DGLAP equations, but the initial condition must be fit to experimental data 
(most importantly from electron-positron annihilations) to extract the non- 
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perturbative information. In this work the DSS [1 14] fragmentation function 
set is usecQ The fact that the cross sections in perturbative QCD can be 
calculated as a convolution of universal parton distribution and fragmentation 
functions (that include the non-perturbative physics) and the perturbatively 
calculable parton level cross section is known as the collinear factorization 
theorem 11161 IT? ]. 

The kinematics of single inclusive particle production is usually written in 
terms of the transverse momentum px and rapidity y of the produced hadron. 
Denoting the center-of-mass energy by a/s, we can define the kinematic 
variables 


x hl 


zP T „y 


Vs 


e , x h2 = 




V~s 


(5.1) 


where Xh 1 is the fraction of the longitudinal momentum of hadron 1 carried 
by the parton taking part in the scattering process originating from hadron 
1, and Xh 2 has the same interpretation for the second hadron. When the 
hadron is produced at central rapidity (y — 0), the situation is symmetric with 
respect to the collided hadrons. When the hadron is produced at forward or 
backward rapidities (at large |y|), we have a large-a; parton from one hadron 
probing the small-a; structure of the second hadron. This discussion is valid 
in two-to-one kinematics where only one parton (which later fragments into 
hadrons) is produced. If higher order processes are considered, one has to 
take into account processes where extra partons are produced in the hard 
scattering process, and in this case X/ ll and Xh 2 become the lower limits for 
the momentum fractions. 


5.2 Particle production cross section from the 

CGC 

The particle production cross section can be computed in two different kine- 
matical limits in the CGC framework. First, we can study the scattering of 
a dilute probe off a dense target and derive the so called hybrid formalism. 
This is applicable in hadron production at forward rapidities, where a large-a: 
parton scatters off a dense target, especially in proton-nucleus collisions where 
“forward” is defined as the proton fragmentation region. The hybrid formalism 
is discussed in more detail in Sec. 15.2.11 

1 Recently it has been pointed out in Ref. j ll5| that the gluon-to-hadron fragmentation 
is too hard in many (NLO) fragmentation functions used in current phenomenological 
applications. 
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Figure 5.1. High-energy quark scattering off the target color field. 


The second kinematical window where we can study the particle production 
is the scattering of two dense objects, such that small -x structure of both 
colliding hadrons is probed. Experimentally this can be realized by measuring 
particle production at central rapidities where the kinematics is symmetric. 
In this case the so called ^-factorization can be used. We will return to this 
in Sec. 15.2.21 


5.2.1 Hybrid formalism 


Let us present a simple derivation for the particle production cross section in 
the CGC framework first derived in Ref. m- 

We consider a quarl0 travelling along the positive light cone with mo¬ 
mentum p. The quark scatters off the target color field, and the final state 
quark momentum is k. The situation is shown schematically in Fig. |5.f| As 


discussed in Sec. 2.2 the incoming quark state can be written as 


|in) = A/j p, i,a) = A f 


d 2 xi 

W) 


■e ip T' XT \p + , x T , i, a). 


(5.2) 


where A f is the normalization factor given in Eq. (2.8) 


When the quark propagates through the target it acquires a phase given 


by the Wilson line, see discussion in Sec. 2.3 The outgoing state reads 


ou 


.t) = A f 


d 2 xt 
(2tt) 2 


e vr "*V ij (x T )\p + ,xr,j, a) 


(5.3) 


The fact that the Wilson line V also depends on the color filed A of the target 
is left implicit. 

To obtain the cross section we calculate the number of quarks in the 
outgoing state by evaluating the expectation value of the quark number 

2 For simplicity we only consider quark scattering here, the derivation for the incoming 
gluon is analogous. 
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operator N q [k ) = b^(k)bjp(k), which measures the amount of quarks with 
momenta k, spin (3 and color j. Summation over repeated indices is left 
implicit. We get 


/ q 2 J2 / 

(2t)2 

x (v} t (x T ')V Ij(iT)) (p + ,i?,j',a\b\j)(.k)bij(k)\p + ,XT,j,a), (5.4) 

where we have taken the average over the possible target color held configura¬ 
tions, denoted by (). 

The number of quarks can be computed by evaluating 


bi,p(k)\p + ,x T , j,a) = 


d 2 yr 


e iyT ' kT bip(k + , y T )\p + i x T ,j, a) 


J (2vr)2 

j^ elVT ' kT biA k+ , VT)bl a {p + , xt )| 0 ) 


= _ e -ixr-hT 5 l Jp,J(p+-k + ), (5.5) 
where the last equality was obtained by applying the anticommutator rule 


(2.6). This result allows us to write 


dAU = V 2 y 

x e - ikT<XT -x T ') [s{p+ _ k + )] 2 5 hj ,5 a p5 l , j 8 a ,p (5.6) 


Recall that px is the transverse momentum of the quark before scattering, 
and kx is the transverse momentum of the scattered quark. The transverse 
momentum of the incoming quark can be taken to be very small compared to 
the transverse momentum transfer in the process, and we can set px = 0. 

The cross section can be computed from the rate of scattering events 
(number of quarks in the final state) normalized by the flux of incoming 
particles (N in = 1 quark in timestep d t, transverse area Sx) and by the 
amount of target particles (here 1): 


d a qA ^ qA 



d N, 


out 


d t 


N • : 

1 v m 

S T dt 


(5.7) 


where we have calculated average over the initial state color and spin, and 
A refers to a proton or a nucleus. The summation over the repeated indices 
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when calculating diV out takes into account the possible spins and colors of the 
final state quark. Calculating the sums we get 


&° qA ^ q+X _ SrAf 2 l§( + _ k+)] 2 [ d 2 x t d 2 x T ' ikT . {xT _ XT ,) 
d?k TM [d{p k J (27r)2(27r)2 e 

x^-(TrV^x T ')V(x T )). (5.8) 

**C 

To evaluate this we change integration variables to r? = Xt — Xt' and 
b t = \{xt + Xt')- The Wilson line correlator (TrEW) now in principle 
depends on three quantities (as one angle can be integrated over due to 
rational symmetry): t't |, |&t| and rx ■ bx- Next, we neglect the rx • bx 
dependence as it should not contribute after we integrate over the impact 
parameter bx- This allows us to write 


da ^g+x ^ S T 6(p + — k + )5(p + = 0) y 
d 3 k (27t) 4 


d 2 b T S(k T ,b T ). 


(5.9) 


The Fourier transform of the dipole scattering matrix is defined as 


S(k Tl b T ) — / dV T e"“^— (WMVyO)) 


(5.10) 


Here the impact parameter dependence is hidden in the target average, as the 
color field density (or saturation scale) can be different at different impact 
parameters. Finally we note that the second delta function gives the size of 
the x~ box, as 5{p + = 0) = f ^e“ '° = Substituting the normalization 


factor from Eq. (2.8) we finally get 


¥k 


(2jt) 


d % T S(k T ,b T ) 


(5.11) 


Let us finally obtain the single particle production cross section for proton- 
proton collisions. As discussed in Sec. 3.3, we assume that the impact 


parameter profile of the scattering matrix S can be factorized, and replace 
f d 2 bx by a 0 /2 which is the proton transverse area measured in deep inelastic 
scattering. The produced particle is assumed to have a large (positive or 
negative) rapidity, and thus one of the incoming protons can be considered 
to be dilute such that the collinear factorization is applicable. As shown in 
Ref. f 1 T9! . the DGLAP-evolution of the parton distribution function /, can 
be included in the above calculation. When calculating the cross section we 
must sum over different parton species. The parton yield is then convoluted 
with the fragmentation function Di^h{z,Q 2 ) which describes the formation 
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of hadron h carrying momentum fraction z of the momentum of the parton i. 
The cross section then becomes 


d (jPP^ h + x 

d 2 k T dy 



~^Xifi(x,Q 2 )Si (^,u) D,^ h (z,Q 2 ). 


(5.12) 


Finally the invariant yield is obtained from the cross section as 


d ]\fpp-+ h + x x <\a pp ^ h+x 
d 2 fc T dj/ (Ti n ei & 2 k T Ay 


(5.13) 


where crj ne i is the total inelastic proton-proton cross section. Here X 2 refers to 
the Bjorken-x at which the target average is performed in Eq. (5.10). The 
dipole amplitude is evaluated in adjoint representation when the gluon channel 
contribution is computed. Note that the factor ~ 0.2 ... 0.3 arises from 

^inel 

the fact that the proton area measured in DIS, cro/2 is different than the 
inelastic proton-proton cross section which relates the invariant yield and the 
cross section to each other. 


When the different areas are properly included, the K factor that is 
sometimes used to scale the LO single inclusive spectra to have the same 
normalization as the data really quantifies how much the leading order 
calculation is off from the data. The yield in proton-nucleus collisions is 
explicitly calculated at a fixed impact parameter and no similar area factor 
appears. 

The hybrid formalism has been successfully applied to phenomenology. 
Combined with the BK evolution, a good description of the RHIC data was 
obtained in Ref. HZDJ. The description was improved in paper [Vj by including 
a proper treatment of different geometric areas in the calculation (the proton 
area measured in DIS and the inelastic proton-proton cross section) to also 
obtain the LO CGC result for the absolute normalization of the spectra. 


This calculation is done at leading order. The next to leading order cross 
section in the hybrid formalism is calculated in Refs. [M , 55J. In principle, 
combining the NLO cross section with the NLO BK evolved dipole amplitude 
(see Sec. 2.5) and NLO parton distribution and fragmentation function allows 
one to calculate the single inclusive cross sections consistently at NLO accuracy. 
However, in Ref. |55j it was shown (without using the NLO evolution of the 
dipole) that the NLO cross section becomes negative. Currently it is not 
known how the NLO calculation should be improved in order to get physically 
meaningful results. 
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5.2.2 Transverse momentum factorization 


Another approach to the gluon production cross section calculation is to use 
so called transverse momentum factorization (hr factorization), which was 
originally derived in Ref. cm and later improved in Ref. [ 122 ] . It can also 
be derived from the Classical Yang-Mills theory, as shown in Ref. [123] , and in 
Ref. [ 124] it was demonstrated numerically that the ^-factorization formula 
agrees with the gluon spectrum computed by solving the Classical Yang-Mills 
equations of motion. For analytical gluon production cross section results in 
this picture we refer the reader to Ref. [125J. 

The ^-factorization formula for the gluon production cross section is 


dcr pp^s+ x 

dyd 2 kTd 2 b T 

where 


2 a s 
Cik 2 T 


d q T d s T 


2 ip(q T , St) <p{k T - qr, b T ~ s T ) 


Qt 


(k T - qr ) 2 


(5.14) 


,2.. . N _ pf ^0/2,4 

“ rvn 


d ~ST<p{k T , st) = 


d 2 r T e 


—ikT'TT 


S(r T ) 


87 r 3 a s ’ VT 

is proportional to the dipole unintegrated gluon distribution: 95 


(5.15) 


k^ iiro '’(k T ), see Eq. ( ]2.28[ ). The dipole amplitude is evaluated in the adjoint 
representation, which is obtained as AAijomtC'C'r) = 2 N(r T ) — N(r T ) 2 ■ To get 
the invariant yield, the cross section is calculated by integrating over 6 ^ and 
the result is divided by the total inelastic cross section. The result is 


d Npp^ +x _ (a 0 / 2) 2 C f 

dyd 2 kT 


d 2 q T 2 


CTinel 87 T 4 k^a s J (27t) : 


q T S(qT)(k T -qT) S(k T ~qT), (5.16) 


where S is the two dimensional Fourier transform of the dipole amplitude 
1 — Na in the adjoint representation. For the proton DIS area we use the 
value cro/2 obtained by fitting the HERA DIS data, see Chapter [3] To get 
the hadron spectra from the parton level yield Eq. ( |5.16[ ) must again be 
convoluted with a fragmentation function. 

Assuming then that the transverse momentum of the produced gluon, |&r|, 


is much larger than the saturation scale the integral in Eq. (5.16) factorizes 


and the hybrid formalism result (5.12) is obtained: 


dN PP^a+x cjq/2 1 


dyd 2 k T CTinel (27 t) 2 
Here the gluon distribution function is 


xg(x, k^)S(kT)- 


(5.17) 


xg(x,k%)= I ^fip(q T ). 


1 0 


Qt 


(5.18) 
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5.3 Proton-nucleus collisions 


To describe proton-nucleus collisions the dipole-nucleus scattering amplitude, 
preferably as a function of the impact parameter, is needed. Here we use 
the dipole amplitude obtained from the fits to HERA electron-proton data 
and generalized to nuclei using standard nuclear geometry as discussed in 
Sec. IH For completeness we write again the dipole-nucleus amplitude at 
x = Xq (which is the initial condition for the BK evolution of the nucleus) 


N A (r T: b T ) = 1 — exp 


Arp ,, To i r TQ 2 s,oV , 

-AT A [b T ) — ---In 


l r r|AQCD 


+ e c 


( 5 - 19 ) 

When the single inclusive yield at fixed impact parameter is calculated 
using the hybrid formalism, one does not get the different area factors in 
Eq. (5.12), and the single inclusive yield is 


dlV(& T ) 

dyd 2 k T 


1 

(2t r) 2 


xg(x,kp)S A (k T ,b T ), 


(5.20) 


where S A is again the Fourier transform of the dipole-nucleus scattering 
matrix in adjoint representation. Similarly when using kx ~factorization only 
one proton area is obtained from the unintegrated gluon distribution of the 


proton, and Eq. (5.16) becomes 


dN pA^g+X Cf(Jo / 2 /• d 2 gT 


dyd 2 k T 


87 i A kpOt. s 


(2w) 


qpS(qT)(k T — q T ) 2 S{kT — qr)- (5-21) 


Physically the reason for having different prefactors in proton-proton and 
proton-nucleus results is due to the consistent usage and separation of the 
two different proton areas: the small- x gluonic area a o/2 determined in DIS 
and the “soft” area of the proton given by the inelastic proton-proton cross 
section cr inel , which gets a large contribution from non-perturbative physics. 

The nuclear effects in proton-nucleus collisions are quantified in terms of 
the nuclear suppression factor R p a, defined as the ratio of particle production 
yields in pA and pp collisions normalized by the number of binary nucleon- 
nucleon collisions Nun- 

dN pA 

^ = < 5 ' 22 > 

The number of binary collisions can be obtained from the nuclear density 
distribution T A using the Optical Glauber model, see Ref. [[HO]- The nuclear 
suppression factor can be interpreted such that if the nuclear effects are 
absent, and the pA collision is the same as a proton colliding with individual 















nucleons, then R p a = 1. Having R p a < 1 can be seen as a signal of saturation 
effects in the nucleus. 

In the high-px limit the nucleus is probed at larger x and at smaller scale 
where the nuclear effects are expected to vanish (except perhaps a small 
difference originating from the isospin symmetry when computing charged 
particle spectra). Let us show that in our framework we naturally obtain 
a value for R ri/ \ which goes to unity at high transverse momenta. First, we 
notice that in this limit the dipole amplitude is evaluated close to the initial 


condition and the single inclusive yield (5.20) becomes 


d N pA 


xg / d 2 r T e 


Hzt'I't 


exp 


AT a 


a 0 N? 


xgATj^NT 


Similarly the proton-proton yield behaves as 

dN pp ~ °^xgN p . 

^"inel 


(5.23) 


(5.24) 


The number of binary collisions is Win = AT/chneb an d thus we find that 
i? pA —>■ 1 independently of the center of mass energy a/s. The fact that 
R p a —t 1 thus originates from a proper treatment of the nuclear and nucleon 
geometry, and it is not a result of the high energy evolution. This can be seen 
from earlier CGC calculations [120, 426) where a significantly larger nuclear 
suppression is predicted for pA collisions, and the main difference between 
this work and Refs. ra and ra is the different treatment of the nuclear 
geometry. 

When calculating the minimum bias particle production yields in proton- 
nucleus collisions we integrate over all impact parameters. At large impact 


parameters the nucleus becomes very dilute, as shown in Sec. 3A, and the 
BIv evolution which is run separately for all impact parameters reduces to 
the BFKL equation and causes the nucleus to grow very rapidly in its edges. 
This unphysical behavior would be cured by an impact-parameter dependent 


BK equation which is not available, see discussion in Sec. 2.4.1 


In practical calculations we use the dipole-nucleus amplitude in the region 
where the saturation scale in the nucleus is larger than that of the proton. 
When the saturation scale falls below that at large | bx\, we use the expanded 


result (5.23) and calculate the yield as 


d N pA = N hin dN pp , 


(5.25) 


which is equivalent to requiring that R p a = 1 at large impact parameters. 

The proton-nucleus collisions can not be divided experimentally into 
the centrality classes similarly as the heavy ion collisions due to smaller 
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p + p—>7 t°/A + X, yfs = 200GeV, K = 2.5 



Figure 5.2. Invariant 7r° and negative 
hadron yields at forward rapidity as a 
function of pt compared with the RHIC 
data |128til30j . Figure from paper (V] , 


p + p —> n°, h* 1 + X, y's = 7000 GeV, K = 1 



Figure 5.3. Invariant charged hadron 
yields at midrapidity as a function of pt 
compared with the LHC data [ 131 L 132] . 
Figure from paper [V]. 


multiplicity and biases introduced by different centrality measures, see e.g. 
Ref. nzzi. The optical Glauber model used in these calculations can be 
used to calculate centrality dependent quantities, but the results can not be 
compared with experimental measurements divided into “centrality classes” 
based on event multiplicity. However, the minimum bias results obtained by 
integrating over the impact parameters can be compared with measurements. 


5.4 Comparison with experimental data 

We calculate in paper [V] single inclusive particle spectra in proton-proton 
collisions using the three different initial conditions for the BK evolution 
fitted to the HERA data in Chapter [3} First, we calculate the neutral pion 
and charged hadron yields using the hybrid formalism and compare with the 


The px dependence of the data is quite well described by the calculation, but 
the absolute normalization must be modified by introducing a factor K = 2.5 
to normalize the calculation to the same level with the data. As the absolute 
normalization is now consistently calculated from the CGC, the K factor 
quantifies how much the leading order calculation differs from the data. 


where the spectra are calculated at midrapidity using /cj- factorization. The 
Pt slope is well reproduced, and no K factor is needed. The unmodified MV 
model initial condition, which gives a relatively good fit to HERA data and 
works with the RHIC measurements, is now clearly not favored by the LHC 


A similar comparison with the LHC data (1311 1132] is shown in Fig. 5.3 


RHIC data |12811130] at a/s = 200 GeV. The results are shown in Fig. 5.2 
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p+Pb/p + p^>h ± + X, y /s = 5020 GeV 



Figure 5.4. Minimum bias nuclear 
modification factor R p a(u = 0) at dif¬ 
ferent centrality classes computed using 
kx factorization and MV 7 , MV e and 
MV model initial conditions compared 
with the minimum bias ALICE data 
[ 133] at smallest px region. Figure from 
paper m- 


p+Pb/p + p-*h* + X,y/i = 5020 GeV 



Figure 5.5. Centrality dependence of 
the midrapidity R p a computed using kx 
factorization and MV e initial condition 
compared with the ALICE data |133] , 
Figure from paper (Vj. 


data. The MV 7 and MV e parametrizations give in practice indistinguishable 
results for the single inclusive yield. 

Let us then study proton-nucleus collisions. First we compute, using 
fcj'-factorization, the nuclear suppression factor R p/ \ in central proton-lead 
collisions at y/s = 5020 GeV and compared with the ALICE minimum bias 
data P331 The results shown in Fig. |5.4| are consistent with the ALICE 
data especially at low px, and all the studied initial conditions for the dipole 
amplitude give the same result for R p a even though the results for the px 
spectra were very different with the MV model initial condition. As the 
systematical uncertainty in the data is relatively large, it can not be said 
whether or not R p a approaches unity in the experimental data. 

The question whether R p a goes to unity at large transverse momenta 
is currently an open question. CMS has observed [134] an enhancement of 
charged particle R p a at transverse momenta > 50 GeV, where R p/ \ ~ 1.4, 
and the result is inconsistent with unity given the experimental uncertainties. 
This behavior is not seen in the ALICE data [ 135] which extends up to 
px ~ 50 GeV. 

Let us then precent our results for the centrality dependence of the nuclear 
suppression factor. First we emphasize again that the centrality classes 
defined using the optical Glauber model do not correspond to the experimental 
centrality classes in pA collisions, and the comparison with the data should 
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p + Pb / p + p —> 7r° + X 



Figure 5.6. Center of mass energy 
(-/s) dependence of the nuclear modifi¬ 
cation factor in neutral pion production 
in minimum-bias p+Pb collisions com¬ 
puted using the MV e initial condition. 
The results at midrapidity y = 0 are 
computed using kx ~factorisation and 
at y = 4 the hybrid formalism is used. 
Figure from paper m 


R vA (j>t = 3 GeV), y/s = 5020 GeV 



Figure 5.7. Rapidity and centrality de¬ 
pendence of the nuclear modification fac¬ 
tor in neutral pion production in 0 — 20% 
most central (solid lines) and minimum 
bias collisions using MV, MV 7 and MV e 
initial conditions. Thin lines at y < 1 
are computed using fcj’-factorization and 
thick lines at y > 2 using the hybrid 
formalism. Figure from paper [Vj. 


be done with extreme care. The midrapidity R p a calculated using the MV e 
initial condition is shown in Fig. 5.5 We find a weak centrality dependence, 
and the results start to deviate only at the most peripheral centrality classes. 
Note that as discussed in Sec. 5.3 we in practice set R p a = 1 in the peripheral 
collisions when the nuclear saturation scale becomes smaller than that of the 
proton, which happens at centralities > 70%. 

To study the center-of-mass energy a/s dependence of the nuclear mod¬ 


ification factor we show in Fig. 5.6 the midrapidity R p a calculated using 
kx~ factorization and R p a(v = 4) calculated using the hybrid formalism at 
two different center of mass energies, corresponding the energies available 
at the LHC before and after the first long shutdown. We find a weak y/s 
dependence. Note that we get explicitly R p a —>• 1 at high p T at all center of 
mass energies. 

Finally we plot in Fig. 5T the rapidity dependence of the nuclear suppres¬ 
sion factor for neutral pions, namely R p a{pt = 3 GeV). As pt is fixed and the 
Bjorken-x probed in the scattering goes like e ~ y , this quantity measures the x 
evolution of the gluon density in the nucleus (where the saturation effects are 
expected to be larger) compared with the proton. The nuclear suppression 
factor close to midrapidity is computed using kx- factorization, and at forward 
rapidities the hybrid formalism is used. The hybrid formalism calculation is 


72 












done using the CTEQ parton distribution function instead of an unintegrated 
gluon distribution function (see discussion in Sec. 2.3), and the quark channel 
is included in the calculation, thus the curve is not exactly continuous. The 
evolution speed of R p a close to midrapidity is slower than at more forward 
rapidities. We also find that all dipole model parametrizations for the initial 
condition give roughly the same evolution speed, and it can be concluded 
that R p a is not sensitive to the details of the initial dipole amplitude. 

We can conclude that R p a does not distinguish between the dipole model 
initial condition parametrizations, as it is mostly sensitive to the BK evolution 
of the dipole amplitude. It would be interesting to compare the CGC calcu¬ 
lations for R p a at different rapidities to experimental data and calculations 
done using different approaches, for example perturbative QCD with nuclear 
parton distribution functions [13611137] . Unfortunately, currently the only 
forward rapidity R p a data comes from RHIC where the lower center-of-mass 
energy limits the utility of the data for this purpose. 

The pt spectra seem to be sensitive to the details of the initial condition 
by clearly disfavoring the MV model. Even though the MV model did 
not give equally good fit to the HERA data as the modified MV model 
parametrizations, comparison with the LHC spectra in proton-proton collisions 
gives an additional constraint which calls for a modification to the MV model 
initial condition. 

In addition to charged particles and pions studied here, it is also possible 
to study for example inclusive J/4/ production cross section. This is done by 
calculating the cross section for cc quark pair production [ 1381 139J , and then 
modelling the cc —> J/T transition using for example nonrelativistic QCD 
(as done e.g. in Ref. [ 140| ) or the so called color evaporation model, where it 
is assumed that fixed fraction of the cc pairs produces the vector meson, and 
the additional color is “evaporated” away as a soft gluon. We have shown in 
Ref. [ 141 j that if the nuclear geometry is properly taken into accoung similarly 
as described in Sec. 5.3 the CGC calculation is in agreement with the LHC 


data in pp and pA collisions. 
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Chapter 6 

Double inclusive scattering and 
correlations 


6.1 Azimuthal correlations as a signature of sat¬ 
uration 


So far in this thesis we have discussed how it is possible to describe single 
particle production in inclusive and diffractive deep inelastic scattering, ul- 
traperipheral heavy ion collisions and hadronic collisions within the CGC 
framework. Going beyond the single inclusive spectrum can provide additional 
tests for the theory, as the multiparton correlations are expected to be more 
sensitive to the detailed dynamics of the colliding particles. 

The azimuthal angle correlations between two pions in the forward rapidity 
have been measured by STAR [ 142 1 1143] and PHENIX [ 129] collaborations at 


RHIC. The striking result from these measurements is that in proton-proton 
collisions the produced hadrons are clearly back-to-back on the transverse 
plane, as expected from naive momentum conservation, but this correla- 



Figure 6.1. Evolution of the parton densities from forward-midrapidity di¬ 
hadron production to the case where both hadrons are produced at forward 
rapidity. 
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Figure 6.2. Dihadron correlation as a function of azimuthal angle difference. 
The trigger and associate particle transverse momenta are shown at the top of 
the panels. Deuteron-gold results are shown for most central (0-20%) and periph¬ 
eral (60-80%) collisions based on charge deposited in the gold-going direction. 
Figure reprinted with permission from Ref. H29|. Copyright (2011) by the Ameri¬ 
can Physical Society. 


tion vanishes when the same measurement is done for central deuteron-gold 
collisions. 

In the kinematical region where one of the pions is produced at midrapidity 
and the second one at forward rapidity (“forward” defined as the direction of 
the projectile), the back-to-back correlation does not disappear. To under¬ 
stand this difference, let us consider the two different kinematical windows 
shown schematically in Fig. 6.1 If one of the produced pions is produced at 


midrapidity, a large amount of longitudinal momentum is needed from the 
target which means that the large-a; structure of the target is probed. On 
the other hand, if both pions are produced at forward rapidity only a tiny 
amount of longitudinal momentum can be taken from the target, and thus the 
significantly more dense small-x part of the target wave function is probed. 
More precisely, if two partons with transverse momenta px,i are produced at 
rapidities y\ and y 2 , then the Bjorken-x of the targed probed in the collision 


is 


x A = 


Pt, ie 


-y 1 




+ 


Pr^e 


~V2 


V~s 


(6.1) 


In the Color Glass Condensate the disappearance of the back-to-back 
correlation peak can be naturally explained. The incoming (valence) quark 
from the projectile emits a gluon either before or after it interacts with the 
target. After the emission the two partons would be back-to-back in the 
transverse plane, but when interacting with the target the partons obtain a 
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Figure 6.3. Quark-gluon production cross section in the scattering of a quark 
off the target. The wavy lines represent multiple interactions with the target. 


transverse momentum kick of the order of the saturation scale Q s . When the 
saturation scale becomes of the same order as the transverse momentum scale 
of the produced partons, the back-to-back correlation is washed out. As the 
saturation scale at the center of the nucleus is much larger than that of the 
proton (or in minimum bias collisions of the nucleus), it is natural to expect 
the correlation to vanish only in central deut.eron-gold collisions. Similarly, 
in this picture we would expect the back-to-back peak to return when the 
transverse momentum of the hadrons is required to be large enough, which is 


exactly what can be seen from the RHIC data in Fig. 6.2 


6.2 Two-particle production cross section 


The cross section for the production of a quark with momentum q and a 
gluon with momentum k, can be computed by calculating the gluon emission 


from a quark as discussed in Sec. 2J2_, and calculating the quark-target and 
gluon-target scattering as in Sec. 2.3 The situation is shown schematically 


in Fig. 6.3, where we note that the gluon can be emitted either before or 
after the interaction with the target. The cross section in coordinate space is 
derived in the CGC framework in Ref. m (see also Ref. [ 1T5] for an original 
CGC derivation in the momentum space): 


foqA^qgX + r dV d 2 X T ' d % T d% 

dk + d 2 k T dq + d 2 q T f ^ ^ J (2n) 2 (2n) 2 (2i r) 2 (27 t) 2 

x e ik *~ b' T )<l>^(x T - bT){S m (bT,x T ,b' T ,x T ') 

ol(3 A 

— S^ipT-, x Ti zt) — S i ' 3 \zt, xt , b' T ) + S(zt, Zt')}- ( 6 . 2 ) 

Here Zt = zxt + (1 — z)br and zt = zxt + (1 — z)b' T and z is the longitudinal 
momentum fraction of the original quark carried by the gluon, z = k + /p + . 
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When calculating the hadron level cross section, the result must be convoluted 
with a parton distribution function and a fragmentation function. In this 
work we use the leading order CTEQ parton distribution function dn and 
DSS fragmentation function [ 114] sets. 


The Wilson line correlators in Eq. (6.2) are 


s^\b T ,b' T ,x T ,x T ') = W. (tv [c(x T )(yt(V)] c ‘‘> 

(6.3) 

S^(b T ,x T ,zr J ) = (Tr (V'(z T ')t‘V(b T )t d ) U*(x T )) (6.4) 

S(z T , V ) = T (Tr (V(z t )V\zt'))) . (6.5) 


Here U again stands for the Wilson line in the adjoint representation and 
( t ) 'ap^ x T ~ M is the wave function for the emission of a gluon with polarization 
A from a quark whose spin changes from a to /9. Note that, as opposite to 
single inclusive cross section and DIS, the cross section can not be expressed 
in terms of the dipole operators (correlator of two Wilson lines) only. The 
different terms in equation (6.2) have a clear physical interpretation: if the 
gluon is emitted before the interaction with the target, both the quark and 
the gluon interact with the target at coordinates Xt and b? (and in the 
complex conjugate amplitude at Xt' and b' T ), and this scattering is described 
by function S^\ If the gluon is emitted after the interaction, only the quark 
is scattering at transverse coordinate zt [zt in the conjugate). The operators 
S<,3) correspond to interference between the two channels. 


In addition to the quark-initiated dihadron production channel, there is 
also a possibility that the incoming parton from the probe is a gluon which 
splits to two gluons or quark-antiquark pair. The cross sections for these 
channels can be calculated similarly and are available in the literature [31., 146] . 
However, as we are here considering only production of two forward rapidity 
hadrons at RHIC, the Bjorken-x for the probe is very large ~ 0.5 and the 
valence quark distribution dominates. If the same analysis is done in case 
of the LHC, the gluon channel should be taken into account. Currently, no 
data on dihadron correlation at forward rapidity at the LHC exists. The 
correlations are measured at midrapidity, but even though the Bjorken-x of 
the target is there similar as in the forward RHIC kinematics, there are other 
diagrams in addition to the ones discussed above that are important in the 
LHC kinematics, see e.g. Ref. (HU- 
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6.3 Double parton scattering 


The expression for the dihadron production cross section, Eq. (6.2), is infrared 


divergent due to a double parton scattering contribution (DPS) included in 
the calculation. In this picture, the DPS contribution corresponds to the case 
where the incoming quark emits a gluon with small transverse momentum 
long before the interaction takes place, and the two parton subsequently 
scatter independently off the target. As the initial transverse momenta of the 
partons are small, the DPS scattering cross section does not depend on the 
azimuthal angle difference A ip between the particles and contributes to the 
two-particle correlation measurement a A ip independent background. 

To see that this contribution is indeed included in the calculation, let us 
consider the case where the quark and the gluon are far away from each other: 
\ut\ = |br — xt\ 3> 1 /Q s , K| = \b'r ~ x t'\ 3> 1 /Q s - Note that the quark 
coordinate in the amplitude and in the complex conjugate can not be far 
away in order to have a non-zero cross section (and similarly for the gluon), 
thus |xt — xt | ~ \br — b ' T | ~ 1 /Q s . We call this case the DPS limit. Now the 


quark and the gluon are not correlated, and the expectation value in Eq. (6.3) 


factorizes. Using the fact that the expectation values must be color singlets 
we obtain 

S (4 \b T , x T , b' T , xt) D ~ g S^ s (b T , x T , b' T , x T ') 


N} - 1 


Tr {V{b T )V\b' T )t d t c )){[U{x T )U\x T ')} cd ) 


N 2 / 1 \ 

^■S(b Tl Hr) [S(;X T , x T ') 2 - W2 j • ( 6 . 6 ) 


N 2 

C 

Here we used the Fierz identities to express the adjoint representation Wilson 
lines in terms of the fundamental representation ones, see discussion in 


Sec. |6.4| This can be interpreted such that the quark scattering is given by 
the correlators S{bT,b ' T ), and 

N? 


N? 




(67) 


is the two point function in the adjoint representation normalized to unity 
at xt = Xt'■ Note that S' 1 - 4 ) is finite in the DPS limit, and the other dipole 
operators S^ and S in Eq. (6.2) vanish. 


In the massless limit the wave function product, that describes the gluon 
emission is 

87T 2 Ut ■ u't 


y, ^a/^T^a/^T) 
a/3X 


k + \ut\ 2 \u' t 


' 12 


(l + (l-z) 2 ). 


( 6 . 8 ) 
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Changing the integration variables to Ut and u' T in Eq. (6.2) it becomes clear 
that the integral over these transverse coordinates is infrared divergent, as 
the correlator is finite in the DPS limit, and a logarithmic divergence is 
obtained from the wave function product. 

Physically, this divergence must be regulated by confinement scale physics 
in the projectile wave function. However, it is not calculable in our pertur¬ 
bative framework. Thus we will subtract the DPS contribution from the 


calculated dihadron production cross section Eq. (6.2) and calculate it sep¬ 


arately by introducing a simple model for the double parton distribution 
function. In practice, we replace the correlator S ^ in Eq. (6.2) by 


^subOr, x T , b' T , x T ') = S {i \b T , x T , b' T , x T ') 


8 f\x T — b T | — - -^ 8 f\ xt ~ b' T \ — — -^ Svl s (b T ,XT,b' T ,x T '), 

V a qcd/ V a qcd/ 


(6.9) 


whith S'ppg defined in Eq. (6.6). Note that with this subtraction the integral 
in Eq. (6.2) vanishes in the DPS limit. It was shown explicitly in paper |llj 


that without this subtraction the cross section (6.2) explicitly reduces to 


a product of quark and gluon scattering matrices and a piece that can be 
identified as the double parton distribution function in the DPS limit. 

To calculate the subtracted double parton scattering contribution we first 
notice that in deuteron-gold collisions there are different processes contributing 
to the Acf) independent background. First, we can have a contribution where 
the two partons are taken from the same proton (or neutron), described by 
the double parton distribution function DPDF. In this case, a kinematical 
constraint x\ + X2 < 1, where X\ and X2 are the momentum fractions of the 
partons, must be satisfied. This is implemented by modelling the double 
parton distribution function (DPDF) following Ref. [148 ]: 


Di?(xi,Xj,Q 2 ) = - 




X-; 


1 - Xi 


+ fi 


Xi 


1 - Xj 




( 6 . 10 ) 


Here i and j denote parton species and the single particle PDFs are evaluated 
at scale Q 2 . 

The second DPS contribution involves taking one parton from the neutron 
and the other one from the proton. In this case the same kinematical constraint 
does not apply, and for the deuteron we can write the DPDF as 

D ij\xi, Xj, Q 2 ) = [fi(xi)f?(xj) + fi(xi)f?{xj)\ . (6.11) 

Here f p is the proton PDF and f n is the same for the neutron. When 
proton-nucleus collisions are considered, this contribution is not included. 
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6.4 Evaluating the multi-point correlators 


In order to calculate the cross section ( 6.2[ ) (with the subtracted DPS con¬ 
tribution from Eq. (6.9)) one has to be able to evaluate the higher point 
functions S'® and S ® written in Eqs. (6.4) and (6.3). 

Using the Fierz identities 


\y{ x T ))ij)\y\ x T)\ki - -^rSuSjk + 2u cd {x T )t c il t d kj 

tij^kl = 2 — 2iV 


( 6 . 12 ) 

(6.13) 


one can show that U C(1 {xt) = 2 TV \V'{xT)t c V{xT)t d \- Using this result, the 
adjoint representation dipoles can be expressed in terms of the fundamental 
dipoles and the higher point functions can be written as 


*S ,(4) (br, x T , b' T , x T ') = - j^S(b T ,b' T ) 

+ S(x T ,x T ')^r (Tr {V{b T )V\b' T )V{x T ')V\x T ))) J (6.14) 
id 

S {3) (b T ,x T ,b' T ) = S(b' T ,XT)S(bT,x T ) - ~^S(b' T ,b T ) 1 (6.15) 


Note that here we have assumed that the expectation value of a product 
of traces factorizes into a product of expectation values. The three-point 
function S^ can be calculated using only the dipole operator (2.24) that 


can be obtained by solving the BK equation. On the other hand, for the 
four-point function an additional correlator of four Wilson lines, known as 
the quadrupole 


Q(x Tl y T ,u T ,v T ) = — (Tr V (x T )V\y T )V {u T )V ] {y T )) (6.16) 


is needed. 

To calculate the quadrupole operator (and its energy evolution) one would 
in principle have to solve the JIMWLK evolution equation. As this would be 
numerically a formidable task, we apply the so called Gaussian approximation 
as presented in Ref. 135. which allows one to write any higher point function 
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in terms of the dipole operators alone. It was shown in Ref. [149J that the 
Gaussian approximation for the quadrupole is close to the solution obtained 
by solving the full JIMWLK equation. Using the Gaussian approximation for 
the quadrupole and taking the large- N c limit we obtain 


S {4 \b T ,x Tl b' T ,XT') = 

_ F{bT-,XT,XT,b ' T ) 

F(bT,x T ',XT,b ' T ) 


S(x T ,x T ') 


S(b T ,x T )S(x T ',b' T ) 


[S(br, x T )S(x T ', b' T ) - S(br, b' T )S(x T ', x T )} , 


where the auxiliary function is 


(6.17) 


F(br, XT-, xt ,b' T ) In S(br, xt) — In S(br, b' T ) + In S(xt, b' T ) — In S(xt, xt) 
F(br, xt ', xt , b' T ) In S(bT, xt) — In SipT, b' T ) + In S(xt ', b' T ) — In S(xt', xt) 

(6.18) 

In previous phenomenological calculations [ 144] 1150] before paper [Hj only 
a so called “naive large- N c limit” for the four point function was used, which 
would give 


S^^bx, XT,b' T , xt) = S(xT,bT)S(xT ,b' T )S(xT,XT'), (6.19) 


which is the same as the first term in the Gaussian approximation (6.17). 


Following Ref. EH we refer to this term as “elastic”. The second term in ( ]6.17[ ) 
is called “inelastic”, and it is exactly this term that remains finite in the DPS 


limit making the integral in Eq. (6.2) infrared divergent. Because the inelastic 


part is not included in previous works ITMIT5U] . the infrared-divergent DPS 
contribution was not found previously. The two-particle correlations were 
also calculated in Ref. (EH by taking the so called back-to-back correlation 
limit where the requirement of having large transverse momenta suppresses 
the DPS contribution. 


6.5 Numerical results 


Let us then present our numerical results for the dihadron correlations at 
RHIC from paper [IT]. The DPS-subtracted cross section (6.2), where the 


quadrupole operator is evaluated using the Gaussian approximation discussed 
in the previous Section, can be evaluated once the dipole scattering amplitude 
N(rx,x) is known. In the calculation we use a BK evolved dipole amplitude. 
For simplicity, and because this work was done before paper where we fit 
the dipole amplitude to the HERA data, the calculation does not use the best 
fit results for the initial condition of the BK evolution. Instead, we just use 
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pp 9 = 2 GeV, = 1 GeV, y\ = y 2 = 3.4 



Figure 6.4. The quark-gluon parton 
level azimuthal correlation in dAu col¬ 
lisions at forward RHIC kinematics. 

A A<p-independent pedestal has been 
added to the “elastic” approximation for 
purposes of visualization. Figure from 

paper ra- 


pp 9 = 2 GeV, p9j? s = 1 GeV, yi = y 2 = 3.4 


p+p, Q^q = 0.2 GeV 2 
nj-Ai i n-.. = (1 79 HpV 2 



12 3 4 


Figure 6.5. The quark-gluon parton 
level azimuthal correlation at forward 
RHIC kinematics in proton-proton, 
proton-nucleus and deuteron-nucleus 
collisions normalized by the pedestal 
contribution. Shown is only the large N c 
result. Figure from paper Bli¬ 


the MV model (2.25) with the initial saturation scale Q S:0 = 0.2 GeV 2 for the 
proton, and scale the saturation scale when changing the proton to a nucleus. 

First we show in Fig. 6A the effect of the “inelastic” contribution to the 
dihadron production cross section (without the DPS contribution). In order 
to get an estimate for the N c suppressed corrections we show both large- N c 
and finite- N c results for the cross section. For comparison, the result obtained 
by using only the elastic part of the correlator (see Eq. (6.17)) is shown. 
Including the inelastic contribution increases the away side peak by a factor 
~ 2 (if the Atp independent pedestal is subtracted). Thus, inclusion of both 
elastic and inelastic contributions is essential for a meaningful comparison 
with the experimental data. 

To study the difference between proton-nucleus, deuteron-nucleus and 
proton-proton collisions we show in Fig. 6J) the two-parton production cross 
section with DPS contribution included divided by the minimum of the Atp 
distribution. Now the larger DPS contribution makes the away side peak 
smaller in deuteron-gold collisions compared to the case where a proton is 
the probe. The disappearance of the away side peak when Q s is increased 
(target is changed from a proton to a heavy nucleus) can be clearly seen 
from the figure. This is a solid prediction from the CGC calculation, and the 
observation of the depletion of the peak in experimental data supports the 
concept of having gluons with larger transverse momenta (of the order of Q s ) 
in the heavy nucleus than in the proton. 
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Figure 6.6. The 7T° azimuthal corre¬ 
lation compared to the preliminary 
STAR [ 143 ] result. The initial satura¬ 
tion scales are Q 2 0 = 1.51 GeV 2 (solid 
line) and Q 2 0 = 0.72 GeV 2 (dashed 
line). Figure from paper [[II]. 


Figure 6.7. Integrated yield under the 
away side peak in central dAu collision 
divided by the corresponding yield in pp 
compared to the forward rapidity part 
of the PHENIX data [129j . Figure from 
paper HH 


To compare with the experimental data we show in Fig. R6 the two-particle 
production yield normalized by the number of trigger particles, known as the 
conditional yield. The calculation is done with two different values for the 
initial saturation scale Q 2 0 of the nucleus calculated by scaling the proton 


saturation scale 0.2 GeV 1 2 by the number of binary collisions in proton- 
nucleus collisions. For the minimum bias collisions we obtain Q 2 0 = 0.72 GeV 2 
and for 0 — 20% centrality class the saturation scale is Q 2 0 = 1.51 GeV 2 . The 
agreement with the STAR data H33] is reasonably good taking into account 
the large experimental uncertainties. Note that the pedestal contribution is 
fixed by hand in order to compare the A ip dependence. As the calculation 
seems to slightly underestimate the peak height, it is clear that using only 
the elastic contribution when evaluating S '* 4 ) would make the agreement with 
the data wors43 


The PHENIX collaboration has measured the nuclear modification factor 
to the integrated conditional yield under the away side peak. We calculate 
this by integrating the area under the peak in dAu collisions and divide the 
result by the same area in proton-proton scattering. The ratio of areas under 


the away side peak, JdAu, is compared with the PHENIX data [129] in Fig. 6.7 


1 A good description with the STAR data by using only the elastic contribution was 

obtained in Ref. CM3 • This agreement was due to a numerical error in the calculation. 
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Data 

Pt range 

pedestal 

exp 

PHENIX pp 

1.1 GeV < Pt 1S < 1.6 GeV 

0.04 

0.095 

PHENIX pp 

1.6 GeV < p T ng < 2.0 GeV 

0.02 

0.079 

PHENIX dAu 

1.1 GeV < Pt 1S < 1.6 GeV 

0.10 

0.176 

PHENIX dAu 

1.6 GeV < Pt S < 2.0 GeV 

0.08 

0.163 

STAR dAu 

2 GeV < p T ng , 1 GeV < pf s < p T ng 

0.02 

0.0145 


Table 6.1. Calculated estimates for the pedestal height compared to the experi¬ 
mental values. In PHENIX results the units are GeV -1 . The dAu values are for 
central collisions where Q 2 0 = 1.51 GeV 2 . Table from paper du. 


The experimental data is measured as a function of 

(p^ ig )e-^ + {p^ s )e~^) 
X f rag = , 


( 6 . 20 ) 


which, at parton level, has an interpretation as the Bjorken-x of the target in 
the event. We calculate this quantity numerically by calculating the ratio of the 
conditional yields in proton-proton and deuteron-gold collisions with different 
hadron momenta and rapidity, and then average the obtained JdAu values 
in each Xf rag bin. The suppression increases as .x'[ r; . lg decreases, as expected 
from the saturation picture, and the theory calculation is in good agreement 
with the PHENIX data. On the other hand, as shown in paper HU, the 
numerical results slightly underestimate the PHENIX data for the conditional 
yield, especially in case of the proton-proton collisions. When calculating 
the conditional yield ratios uncertainties in e.g. single inclusive baseline 
calculations partially cancel, which makes the JdAu calculation perhaps more 
solid than the results for the conditional yields. 

Finally we calculate an estimate for the Ap independent pedestal contri¬ 
bution by summing the DPS contribution and the Aip independent part of the 


clihadron production cross section (6.2). The results are shown in Table 6.1 


and are compared with the experimental data. As our model to calculate the 
DPS contribution is quite rough and there are theoretical uncertainties also 
related to the calculation of the single inclusive spectra, we do not expect to 
get a perfect match with the data, the results being off by roughly a factor of 
2 . 


To conclude, we emphasize that the disappearance of the away side peak 
in central deuteron-gold collisions compared with the proton-proton ones is 
naturally described within the CGC framework both qualitatively and (at least 
semi) quantitatively. In the CGC calculations the inclusion of the “inelastic 
contribution” is essential, and it is also necessary to correctly separate the 
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double parton scattering contribution. 

A measurement of the correlation of two forward rapidity particles at the 
LHC would provide interesting additional set of data. The advantage of the 
LHC is that with significantly higher center-of-mass energy the two-particle 
production takes place fare away from the kinematical boundary, in contrast 
to the RHIC kinematics. This would require an inclusion of the gluon channel 
in the calculation, as the projectile proton would also be probed at relatively 
small x. This channel is calculated in Refs. \m ma, but no numerical 
predictions exist so far. 

In addition to correlations of two hadrons, a powerful probe of QCD 
dynamics could be photon-hadron correlations at forward rapidities. The 
photon could be used to fix the parton level kinematics, in contrast to dihadron 
production where the hadronization procedure smears the process. The second 
advantage would be that no higher-point functions are needed to evaluate 
the cross section, which makes the theory calculation more solid as one does 
not have to use e.g. Gaussian approximation to approximate the quadrupole 
operator. For CGC predictions of photon-hadron correlations, the reader is 
referred to Refs. [1521H53] . 






Chapter 7 

Conclusions and outlook 


The work presented this thesis shows, when combined with the other works 
cited in the text, that the Color Glass Condensate is in agreement with the 
currently available experimental data. The picture can also be systematically 
improved in the future. 

Particle production cross sections calculated for inclusive and exclusive 
processes in papers cheiei are compatible with the available experimental 
data. The precise combined HERA deep inelastic scattering measurements 
are in agreement with the CGC calculations; taking the non-perturbative 
input from the DIS measurements, it is possible to obtain a good description 
of single particle production in proton-proton and proton-nucleus collisions 
(paper fV] and Chapter [5]). However, in perturbative QCD calculations the 
next to leading order corrections to the cross section are known to be large 
and so far we have only done calculations at leading order. It is therefore 
an important goal for CGC phenomenology to calculate the single particle 
spectra and nuclear suppression factors at next to leading order accuracy. 
First steps in this direction have been taken, as discussed in Chapter [5j 

In addition to single particle spectra, two-particle correlations can be used 
to obtain more detailed information of the strong color fields of the nucleus at 
small-ax The fact that the RHIC forward dihadron correlation measurements 
are naturally explained qualitatively, and even semi-quantitatively, in the CGC 
framework (as discussed in paper [Hj and in Chapter [6]) strongly suggest that 
we are able to see the gluon saturation already at RHIC energies. Measuring 
the correlation at the LHC as well would be extremely interesting, since at 
the LHC we would be far away from the kinematical boundary, in contrast to 
RHIC. A much cleaner process would be photon-hadron correlations, and if 
an electron-ion collider will be built in the future, it will also open many new 
possibilities for correlation measurements. 

Exclusive vector meson production, as discussed in Chapter [3J can be 
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used as a direct probe of the gluon distribution. If an electron-ion collider 
will be realized, diffractive production of vector mesons can be used to probe 
gluon densities and density fluctuations in the nucleus, giving access to 
spatial distribution of gluons as well. Towards this goal, we have calculated 
cross sections for incoherent diffractive vector meson production in paper fl] 
and studied the centrality selection in paper [VI]. Our proposal for the 
centrality estimator is the multiplicity of the “ballistic” protons produced in 
the scattering process. Inclusion of a more detailed description for the nuclear 
breakup has a potential to make this estimator more accurate. Incoherent 
and coherent diffraction can also be studied with the current experimental 
facilities in ultraperipheral heavy ion collisions, as discussed in Sec. 4.5 In 


paper HV| we emphasized the importance of the simultaneous description 
of both coherent and incoherent vector meson production when comparing 
the theory calculations with the data. An important next step for the CGC 
calculations in this held is to figure out how to do these calculations properly 
with a BK evolved dipole amplitude and simultaneously describe production 
of different mesons. 

An important ingredient in the next to leading order CGC phenomenology 
will be a small-x evolution equation at next to leading order accuracy. Towards 
this goal, we have made a detailed analysis of the NLO BK equation 


see 


Sec. 2.5 and paper EU) and proposed a method to include running coupling 
corrections in the JIMWLK equation (Sec. 2.4.2, paper (III]). The NLO 
BK equation is shown to be unstable when solved with phenomenologically 
relevant initial conditions, which suggests that the equation would require 
a proper resummation of higher order contributions before it is ready for 
applications. 

Signatures of the saturation phenomena described by the Color Glass 
Condensate are seen in the experimental data and the CGC calculations are in 
qualitative (and in many cases quantitative) agreement with the measurements. 
Taking the theory to next to leading order accuracy, while simultaneously 
comparing the saturation physics calculations with new interesting measure¬ 
ments from current and especially future accelerators, is an important and 
compelling task for the community. 
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